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ABSTRACT:  The  Army  maintains  an  aging  inventory  of  over  143,000  structures,  many  of  which  are  in  urgent  need 
of  maintenance  and  repair.  In  some  instances,  these  aging  structures  fail  to  meet  prevailing  seismic  engineering 
codes.  In  the  United  States  alone,  30  percent  of  the  Army’s  structures  use  unreinforced  masonry  (URM)  walls.  As 
such,  URM  structures  have  inadequate  plane  lateral  strength  and  are  prone  to  failure  during  seismic  events.  Upgrad¬ 
ing  these  structures  to  meet  existing  seismic  codes  often  requires  the  use  of  new  materials  and  systems  such  as  fiber- 
reinforced-polymer  (FRP)  composites.  The  versatility  and  resilience  of  such  advanced  composites  make  them  ideal 
candidate  materials  for  reducing  the  cost  of  seismic  rehabilitation  of  DOD  facilities.  Though  many  studies  have 
demonstrated  the  effectiveness  of  FRP  reinforcements  as  structural  upgrades  for  masonry  walls,  little  has  been  done 
to  explain  their  impact  on  the  building  envelope.  This  study  was  undertaken  to  discover  and  define  the  combination 
of  building  envelope  and  hygrothermal  conditions  that  might  result  in  vapor  liquefaction  at  the  ceramic-epoxy  inter¬ 
face,  and  subsequently  debond  or  delaminate  the  applique.  This  research  developed  fundamental  models  and  asso¬ 
ciated  material  parameters  to  predict  thermal  and  moisture  transport  across  dissimilar  building  materials — 
specifically  FRP  composite  appliques  to  concrete  masonry  units  (CMUs) — used  for  seismic  upgrades. 
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Conversion  Factors 


•  m  m  ® 

Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  as 
follows: 


Multiply 

By 

To  Obtain 

Acres 

4,046.873 

square  meters 

cubic  feet 

0.02831685 

cubic  meters 

cubic  inches 

0.00001638706 

cubic  meters 

degrees  (angle) 

0.01745329 

radians 

degrees  Fahrenheit 

(5/9)  x  (°F  -  32) 

degrees  Celsius 

degrees  Fahrenheit 

(5/9)  x(°F- 32) +  273.15. 

kelvins 

Feet 

0.3048 

meters 

gallons  (U.S.  liquid) 

0.003785412 

cubic  meters 

horsepower  (550  ft-lb  force  per  second) 

745.6999 

watts 

Inches 

0.0254 

meters 

kips  per  square  foot 

47.88026 

kilopascals 

kips  per  square  inch 

6.894757 

megapascals 

miles  (U.S.  statute) 

1 .609347 

kilometers 

pounds  (force) 

4.448222 

newtons 

pounds  (force)  per  square  inch 

0.006894757 

megapascals 

pounds  (mass) 

0.4535924 

kilograms 

square  feet 

0.09290304 

square  meters 

square  miles 

2,589,998 

square  meters 

tons  (force) 

8,896.443 

newtons 

tons  (2,000  pounds,  mass) 

907.1847 

kilograms 

Yards 

0.9144 

meters 

Systeme  International  dVnites  (“International  System  of  Measurement”),  commonly  known  as  the  “metric  system.1 
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1  Introduction 


Background 

The  Army  maintains  an  aging  inventory  of  over  143,000  structures,  many  of 
which  are  in  urgent  need  of  maintenance  and  repair.  In  some  instances,  these 
aging  structures  fail  to  meet  prevailing  seismic  engineering  codes.  In  the  United 
States  alone,  30  percent  of  the  Army’s  structures  are  constructed  with  unrein¬ 
forced  masonry  (URM)  walls.  URM  structures  have  an  inadequate  plane  lateral 
strength  and  are  consequently  prone  to  failure  when  subjected  to  natural  haz¬ 
ards  (e.g.,  seismic  events). 

Upgrading  these  structures  to  meet  existing  seismic  codes  often  requires  the  use 
of  new  materials  and  systems  such  as  fiber-reinforced-polymer  (FRP)  composites. 
FRP  and  other  advanced  composite  materials  offer  many  advantages  over  more 
conventional  building  materials.  Their  high  specific  stiffness  and  specific 
strength  properties,  coupled  with  their  intrinsic  corrosion  resistance  and  ease  of 
application  make  them  ideal  for  incorporation  into  many  existing  lightweight 
composite  structural  elements.  They  can  conform  to  non-planar  surfaces,  can  be 
cut  and  applied  in  various  forms,  and  if  applied  properly,  need  not  cover  the  en¬ 
tire  surface  of  the  URM  wall.  When  bonded  to  URM  walls,  they  can  greatly  en¬ 
hance  the  ductility  of  the  wall.  Moreover,  many  of  these  materials  are  “environ¬ 
mentally  friendly.”  They  consume  less  energy  during  fabrication  than  do 
traditional  materials,  and  can  be  applied  in  situ.  These  features  combine  to  re¬ 
duce  construction  and  shipping  costs,  which  becomes  particularly  evident  (and 
desirable)  when  materials  are  required  in  remote  construction  locations.  In 
summary,  the  versatility  and  resilience  of  advanced  composites  make  them  ideal 
candidate  materials  for  reducing  the  cost  of  maintenance  and  for  seismic  reha¬ 
bilitation  of  DOD  facilities. 

Though  many  studies  have  demonstrated  the  effectiveness  of  FRP  reinforce¬ 
ments  as  structural  upgrades  for  masonry  walls  (Al-Chaar  and  Hassan  2002; 
Marshall,  et  al.  1999;  Saadatmanesh  1997;  Kim  and  Lee  1998;  Sakata  1983;  An- 
and  and  Gandhi  1983),  little  has  been  done  to  explain  their  impact  on  the  build¬ 
ing  envelope.  A  crucial  feature  in  the  use  of  any  composite  applique  is  the  ability 
of  the  composite  to  remain  in  intimate  structural  contact  with  the  supporting 
member  undergoing  retrofit.  Detailed  knowledge  of  the  physical  parameters  of 
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both  applique  and  substrate,  compounded  with  any  environmental  factors,  will 
collectively  define  the  parameter  space  of  application.  It  is  of  particular  interest 
to  model  the  moisture  migration  in  an  FRP  upgraded  URM  wall  structure  since 
the  adhesive  interface  between  the  FRP  laminate  and  the  host  masonry  can  be 
susceptible  to  moisture  entrapment.  Moisture  entrapment  can  lead  to  a  reduced 
load  transfer  at  the  FRP/masonry  interface  or  result  in  premature  FRP  debond¬ 
ing.  Thus,  the  mechanical  durability  of  the  composite-substrate  interface  and 
attendant  viscoelastic  sensitivities  to  moisture  and  temperature  fluctuations  can 
be  the  defining  parameters  for  many  applications.  This  matter  is  further  com¬ 
plicated  by  the  fact  that  composites  may  be  exposed  to  many  different  harsh  ex¬ 
ternal  environments  ranging  from  extremely  cold  climates  (e.g.,  Fairbanks,  AK) 
to  warm,  humid  climates  (e.g.,  Key  West,  FL). 

A  decade  of  field  data  and  laboratory  research  has  continued  to  justify  and 
champion  the  use  of  FRP  upgrades  for  URM  when  executed  in  a  consistent  and 
conscientious  manner  by  trained  professionals.  However,  research  has  shown 
that  the  true  kinetics  of  moisture  absorption  (FRP/URM)  under  freeze-thaw  ex¬ 
posure  in  aqueous  solutions  is  not  yet  well  understood.  The  current  body  of  lit¬ 
erature  does  not  explain  the  presence  of  microcrack  and  fiber-matrix  debond  re¬ 
lated  degradation  and  possible  mass  dislocation  under  these  conditions  (Rivera 
and  Karbhari  2002).  As  such,  long-term  predictions  made  using  models  based  on 
short-term  phenomena  of  several  studies  may  lead  to  an  overly  optimistic  de¬ 
ployment-parameter  envelope. 

During  the  same  period,  computer  technology  made  significant  strides.  The  cost 
of  computer  simulations  declined  exponentially  with  the  appearance  of  inexpen¬ 
sive  memory  and  faster  clock  speeds.  User-friendly  software  has  kept  pace  with 
these  rapid  hardware  innovations,  particularly  those  employing  Finite  Element 
Codes.  At  the  same  time,  a  branch  of  research  known  as  “Building  Physics”  was 
begun  in  Northern  Europe  and  Canada.  Careful  experiments  on  common  build¬ 
ing  materials  using  innovative  measurement  techniques  generated  research  data 
that  explained  how  temperature  and  moisture  commingled  in  a  collective  migra¬ 
tion  through  test  specimens. 

The  coincidence  of  these  events  has  made  computer  simulation  a  valuable  re¬ 
search  tool  to  help  resolve  the  role  of  hygrothermal  effects  in  the  application  of 
FRP  composites  to  the  building  envelope.  This  study  was  undertaken  to  discover 
and  define  the  combination  of  building  envelope  and  hygrothermal  conditions 
that  might  result  in  vapor  liquefaction  at  the  ceramic-FRP  interface,  and  subse¬ 
quent  debonding  of  the  applique. 
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Objective 

The  objective  of  this  research  was  to  develop  fundamental  models  and  associated 
material  parameters  to  predict  thermal  and  moisture  transport  across  dissimilar 
building  materials — specifically  FRP  composites  appliques  to  concrete  masonry 
units  (CMUs) — used  for  seismic  upgrades. 


Approach 

This  work  involved: 

1.  Developing  a  detailed  description  of  the  physical  mechanisms  responsible  for 
(non-saturated)  moisture  migration  through  porous  materials,  and  of  the  related 
role  of  temperature  in  this  transport  process. 

1.  Performing  a  detailed  literature  review  to  obtain  the  physical  parameters  that 
appear  in  the  coupled  hygrothermal  equations  (which  have  both  sensitivity  to 
moisture  and  temperature),  and  details  of  how  these  parameters  are  measured. 

2.  Modification  of  commercial  computer  software  to  solve  the  coupled  set  of  mois¬ 
ture  and  thermal  differential  equations  derived  in  the  previous  steps. 

3.  Using  these  derived  results  to  provide  modified  strain  boundary  conditions  as 
initial  data  required  in  standard  stress-strain  analyses  of  mechanical  systems. 

4.  Interpreting  the  hygrothermal  and  mechanical  results  to  provide  conclusions  and 
recommendations  for  the  deployment  of  FRP-CMU  appliques  in  overall  building 
envelope  design. 


Mode  of  Technology  Transfer 

It  is  anticipated  that  the  modeling  effort  developed  in  this  work  will  be  used  in 
conjunction  with  mechanical  strain  models  to  improve  the  placement  and  per¬ 
formance  of  FRP  composites,  and  ultimately,  to  improve  overall  building  enve¬ 
lope  design.  Immediate  plans  are  to  incorporate  the  fundamental  models  and 
results  from  this  basic  research  project  into  seismic  upgrade  AT41  projects,  in¬ 
cluding  “Accelerated  Performance  and  Weathering  of  Structural  Repair  and  Up¬ 
grade  Composites,”  and  “Seismic  Rehabilitation  of  URM  Walls.”  The  results  of 
this  work  will  continue  to  be  disseminated  through  scientific  and  trade  litera¬ 
ture.  Current  journal  articles  and  professional  presentations  resulting  from  this 
work  unit  are: 

Lin,  M.W.,  M.  Kahbahda,  J.  Berman,  C.A.  Feickert,  and  A.O.  Abatan.  2002.  “Nonlinear  Finite 
Element  Modeling  of  Moisture  Migration  in  a  Masonry  Structure  with  FRP  Upgrade,” 
ISBN  0864912323,  presented  at  CanCNSM,  Vancouver,  B.C.,  19-23  June  2002. 
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A.  Abatan,  M.  Lin,  et  al.,  “Durability  of  Fiber  Reinforced  Polymer  Composites  Under  Varying 
Environmental  Behavior,”  The  4th  Canadian  International  Composites  Conference, 
CANCOM  2003,  19-22  August  2003,  Ottawa,  Canada. 

This  report  will  be  made  accessible  through  the  World  Wide  Web  (WWW)  at 
URL: 

http://www.cecer.armv.mil 
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2  Moisture  Infiltration  of  Porous 
Structures 


During  the  past  several  decades,  the  area  of  building  physics  research  has  devel¬ 
oped  into  a  mature  multifaceted  discipline  particularly  in  the  countries  of  North¬ 
ern  Europe.  The  discipline  owes  its  rigorous  heritage  to  the  formal  investiga¬ 
tions  of  earlier  researchers  like  Philip  and  De  Vries  (1957),  and  Luikov  (1975), 
who  applied  thermodynamics  and  Navier-Stokes  considerations  to  the  detailed 
theory  of  heat  and  mass  transfer  in  porous  materials.  Such  results  have  enabled 
engineers  to  perform  detailed  modeling  calculations  on  common  building  materi¬ 
als  and  thereby  to  more  accurately  assess  their  suitability  to  a  given  moisture 
and  temperature  regime. 

Many  of  the  appliques  and  substrates  in  particular  tend  to  be  porous,  and  to 
have  material  properties  that  are  particularly  sensitive  to  ambient  humidity 
(moisture)  and  temperature.  Porous  materials  tend  to  absorb  moisture  from  the 
surrounding  moist  air.  This  uptake  of  water  vapor  (moisture)  is  referred  to  as 
“sorption,”  which  continues  until  a  state  of  equilibrium  is  reached.  The  moisture 
content  of  the  material  at  this  state  is  referred  to  as  the  “equilibrium  moisture 
content,”  characterized  by  the  parameter  6.  The  moisture  content,  9,  is  typically 
measured  in  units  of  mass  per  unit  volume  of  the  dry  material  (kg/m3,  the  unit 
used  in  this  work),  u  the  mass  of  adsorbed  water  per  mass  of  the  dry  material 
(kg/kg),  or  sometimes  volume  of  water  per  volume  of  dry  material  (m3/m3).  The 
moisture  content,  6,  is  related  to  these  other  measures  by  multiplying  them  by 
the  (dry)  density  of  the  porous  material  or  the  density  of  water,  respectively. 

Because  of  this  sorptive  property,  there  exists  an  important  defining  relationship 
between  the  vapor  pressure  of  water,  or  relative  humidity  (RH)  of  the  ambient 
air,  and  the  equilibrium  moisture  content  within  the  material.  The  essence  of 
this  correlation  is  described  by  a  curve  called  the  “sorption  isotherm.”  Figure  1 
shows  a  typical  sorption  isotherm,  characteristic  of  a  large  class  of  porous  build¬ 
ing  materials  (Valen  1998;  van  Zijl  1999;  Janz  1997).  The  curve  relates  the 
moisture  content  6,  to  ambient  humidity  h,  and  has  a  characteristic  sigmoid 
shape,  with  the  slope  of  the  linear  region  defining  the  moisture  capacity  0  of  the 
material.  At  constant  relative  humidity  of  the  surroundings,  the  equilibrium 
moisture  content  decreases  with  increasing  air  temperature,  however  the  differ¬ 
ence  in  moisture  content  is  small  (Valen  1998). 


6 


ERDC/CERL  TR-03-20 


Genaric  Sorption -Desorption  Isotherms 


Figure  1.  Atypical  sorption-desorption  isotherm  for  a  porous  cement  based  material,  defining 
the  hygroscopic  range  and  maximal  moisture  content  beyond  98%  RH. 

The  figure  defines  the  hygroscopic  range,  critical  moisture  content,  capillary  and 
maximum  moisture  content  of  a  porous  material.  The  range  of  RH  between  0 
and  98  percent  is  defined  as  the  hygroscopic  range.  Beyond  the  hygroscopic 
range,  when  the  RH  is  higher  than  98  percent,  capillary  moisture  condensation 
occurs  within  the  pores.  At  this  stage,  there  is  as  yet  no  continuity  of  the  liquid 
in  the  capillaries.  This  process  continues  until  a  critical  moisture  content,  Gcr 
(defined  as  the  lowest  moisture  content  necessary  to  initiate  contiguous  moisture 
transport  in  the  liquid  phase),  is  established.  Below  this  level,  moisture  is 
transported  only  in  the  vapor  phase,  even  though  there  exists  (isolated)  capillary 
moisture  condensation.  In  addition,  because  of  the  sponge-like  feature  of  porous 
materials  (i.e.,  the  porosity  is  fractal),  the  reverse  process  of  moisture  outgasing 
is  characterized  by  hysteresis,  thereby  ensuring  a  nonlinear  diffusion  mecha¬ 
nism.  For  the  purposes  of  this  study: 

1.  All  moisture  (liquid  and  vapor)  is  assumed  to  be  transported  within  the  hygro¬ 
scopic  range.  No  freezing,  boiling,  or  melting  occurs. 

2.  Both  adsorption  and  desorption  isotherms  are  single  valued  functions  of  the  am¬ 
bient  relative  humidity.  Consequently,  the  models  presented  here  only  describe 
monotonic  wetting  or  drying  processes.  Cyclic  excursions  within  the  hygroscopic 
region  require  a  more  detailed  analyses  and  knowledge  of  the  material  proper¬ 
ties. 
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Moisture  and  Thermal  Transport  Equations 

To  quantify  the  migration  process  further,  one  applies  Fick’s  law  for  moisture 
flux  and  the  general  continuity  equations  to  a  given  media.  The  resulting  cou¬ 
pled  equations  assume  that  the  moisture  content,  0,  is  not  lost  through  any  cur¬ 
ing  mechanism,  but  serves  (generally,  along  with  T)  as  the  driving  potential  (i.e., 
state  variables)  for  diffusion.  Finally,  gravitational  and  dynamic  pressure  differ¬ 
ence  (i.e.,  wind)  effects  will  be  considered  negligible  for  the  particular  formula¬ 
tion  given  below. 

f)  f) 

—  +V0^  =  O  Eq  1 

dt  * 

cf  +  DfNO  +  ZJt’VT  =  0 

Here  q  depicts  the  flux  vector  in  units  of  (kg/m2.s)  and  the  functions  Do  and  Dt 
refer  to  the  total  (liquid  plus  vapor)  moisture  diffusivity  and  thermal  (liquid  plus 
vapor)  diffusivity,  respectively.  The  units  of  Do  generally  depend  on  those  of  9, 
choosing  0  in  units  of  (kg/m3)  yields  units  of  (m2/s)  for  Do,  with  corresponding 
similar  considerations  for  Dt. 

The  compactness  of  these  equations  belies  considerable  subtly  in  their  formula¬ 
tion  (Janz  1997),  and  implicitly  incorporates  the  fact  that  present  technology 
cannot  make  a  meaningfully  distinction  between  in  situ  moisture  in  its  liquid  or 
vapor  phases.  During  any  water  transport  process,  there  will  always  exist 
evaporation  or  condensation  at  the  water  menisci.  Thus,  there  will  be  a  phase 
transformation  between  the  vapor  and  liquid  phases.  The  equations  that  de¬ 
scribe  each  of  these  separate  phase  flows  contain  additional  terms  that  sum  to 
zero  when  the  total  moisture  is  considered,  as  in  Equation  1  above. 

Analogous  formulations  of  Equation  1  exist  depending  on  the  details  and  choice 
of  driving  potential  (Valen  1998),  and  are  normally  used  within  the  hygroscopic 
region  (0-98  percent  RH).  However  the  analytic  form  of  Equation  1  is  sufficiently 
general  so  as  to  be  applicable  to  the  whole  region  from  dry  material  up  to  100 
percent  humidity,  provided  modifications  to  Do  and  Dt  are  made  that  incorporate 
the  new  hygro-thermal  physics  that  is  present  when  h  >  98  percent  RH.  Fur¬ 
ther,  Equation  1  is  implicitly  formulated  with  the  notion  that  all  moisture  trans¬ 
port  of  capillary  absorption  measurements  occurs  in  liquid  phase,  and  all  mois¬ 
ture  transport  in  the  cup-method  measurement  occurs  in  vapor  phase.  In 
practice,  it  is  more  appropriate  to  use  dry-cup  values  and  adjust  the  liquid 
transport  coefficients  to  include  surface  diffusion.  A  notable  exception  are  plas¬ 
tic  film  coatings,  where  the  vapor  permeability  itself  can  increase  considerably 
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with  humidity  because  water  molecules  creep  between  the  polymer  chains  and 
thus  change  the  polymer  structure  (Kuenzel  et  al.  2000). 

In  addition,  it  is  also  important  to  keep  in  mind  the  mathematical  model  for 
which  the  measured  data  will  be  used.  For  example,  transport  coefficients 
measured  under  isothermal  conditions  cannot  automatically  be  used  in 
non-isothermal  models.  If  the  moisture  content  9  and  the  temperature  T  are 
used  as  state  variables  in  the  non-isothermal  model,  the  transport  coefficients 
will  be  De(9,  7)  and  Dt{T,  9)  respectively.  In  this  case,  measurements  must  be 
performed  at  varying  moisture  contents  and  moisture  gradients,  and  at  different 
temperature  levels  and  temperature  gradients  (Claesson  1993).  In  doing  so,  one 
is  then  able  to  determine  the  first  order  contributions  of  a  general  Taylor  series 
expansion  for  Dd,9,  7)  and  Dt{T,  9).  Fortunately,  for  many  practical  instances, 
one  may  simply  consider  De  =  Dd.8),  with  the  T  sensitivity  being  a  higher  order 
term.  Similarly  Dt  is  considered  to  have  a  slowly  varying  sensitivity  to  the  local 
temperature  and  relative  humidity. 


Interfacial  Conditions  of  Continuity 

Equations  1  are  incomplete  as  presented;  all  differential  equations  require  a  set 
of  boundary  conditions  for  their  complete  solution.  In  most  cases,  these  condi¬ 
tions  are  defined  by  some  combination  of  relevant  physics  and  constraining  ge¬ 
ometry.  These  constraints  are  typically  developed  and  executed  at  the  surfaces 
or  boundaries  of  the  system,  which  gives  rise  to  the  nomenclature.  Two  distinct 
materials  that  are  in  close  contact  with  each  other,  and  share  a  common  surface 
with  some  degree  of  interpenetration,  are  said  to  be  in  “Natural  Contact.”  As 
such,  the  moisture  content,  9,  of  the  two  materials  can  have  vastly  different  val¬ 
ues,  so  will  be  discontinuous  at  the  interface.  However  for  a  system  that  is  in 
equilibrium  at  any  instant,  the  total  (entering  and  leaving)  moisture  flux  per¬ 
pendicular  to  the  surface  of  a  given  volume  must  remain  constant.  Furthermore, 
it  is  an  experimental  fact  that  the  vapor  pressure  of  water  (related  to  the  relative 
humidity)  is  the  same  at  both  surfaces.  For  such  surfaces  in  natural  contact,  the 
continuity  of  normal  flux  vector  component  (Equation  1),  qn,  and  the  relative 
humidity,  h,  across  an  (A-B)  interface: 

qnA  =  qnB  Eq  2 


hi  =  hs 


provide  the  required  set  of  boundary  conditions  for  Equation  1. 
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Thermal  Transport  Equations 

The  set  of  Equations  1  and  2,  are  insufficient  to  completely  describe  moisture 
migration  through  any  real  material  system  with  a  finite  specific  heat  capacity, 
c,  since  such  moisture  transport  intrinsically  involves  a  two-phase  (thermal)  sys¬ 
tem.  A  second  set  of  equations  that  describe  the  thermal  fluxes  through  the  in¬ 
dividual  system  components  (sometimes  is  referred  to  as  the  “energy  equations”) 
is  required: 

pc{dTldt)  =  VO (JctVT)  +  hvVU(DTVT)  Eq  3 

Here: 

p  =  the  host  material’s  mass  density  in  kg/m3,  and  is  typically  a  function  of 
the  amount  of  moisture  contained  within  a  porous  material. 

c  =  the  host  material  specific  heat  in  J/kg.K,  and  is  also  a  function  of  mois¬ 
ture  content. 

kr  =  the  host  material  thermal  conductivity  in  Watts/m.K,  and  is  similarly  a 
function  of  moisture  content. 

hv  =  the  latent  heat  of  evaporative  phase  change,  in  J/kg.  Handbook  values 
for  water  are  listed  as  2,500.9  kJ/kg  at  0  °C,  and  2,256.8  kJ/kg  at 
100  °C. 

For  heat  and  moisture  transport,  the  traditional  boundary  conditions  are  re¬ 
quired,  in  particular  the  continuity  of  normal  component  of  thermal  flux,  and 
temperature  at  an  interface.  Dt  is  as  defined  in  Equation  1  above,  however  a 
simplification  is  possible.  Explicit  calculation  of  using  the  formalism 

found  in  Appendix  A  clearly  shows  that  this  term  is  several  orders  of  magnitude 
smaller  than  a  typical  kr.  However,  the  term  has  been  included  in  Equation  3 
for  completeness. 

Collectively,  Equations  1  to  3  are  sufficient  to  define  moisture  migration  within 
the  hygroscopic  regions  (Figure  1)  for  each  set  of  materials  in  a  masonry-FRP 
applique.  In  the  most  general  case,  these  nonlinear  equations  must  be  solved 
iteratively  with  T  as  the  connecting  variable.  However  for  the  ceramic-epoxy  - 
FRP  system,  the  kr  of  a  typical  CMU  is  at  least  an  order  of  magnitude  greater 
than  that  of  the  epoxy  or  FRP.  This  affords  a  great  simplicity  by  effectively  de¬ 
coupling  the  two  equations  in  temperature  but  not  necessarily  in  time.  In  doing 
so,  one  is  able  to  first  solve  Equation  3  and  thereby  prescribe  the  thermal  field 
distribution  for  a  given  moisture  profile,  B  =  B(r,  t  =  constant).  The  resulting 
thermal  solution,  T=  T[r,  t-  constant),  can  then  be  used  to  solve  for  a  new  sub¬ 
sequent  moisture  migration  resulting  from  the  prescribed  thermal  field. 
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Some  Particulars  for  D#(0)  in  Cementitious  Materials 

Figure  2  shows  the  general  features  of  a  diffusivity  curve,  Do  =  DdO)  for  several 
common  ceramic  materials.  The  actual  details  depend  on  specific  material  prop¬ 
erties.  However  some  generalizations  are  possible: 

logio(D9  m2/s])  »A9  +  B  (dimensionless). 

for:  0.0  <  9  (kg/m3)  <  250,  and  Eq  4 

0.013  <  A  (m3/kg)  £  0.022  ;  -11.0  <  B  Z  -9.0 

These  values  of  A  and  B  were  summarized  from  van  Zijl  (1999)  and  characterize 
materials  such  as  mortar,  sand-lime  1,  clay  brick,  and  calcium  silicate.  Equation 
4  provides  a  reasonable  description  of  De  for  these  materials  during  the  sorption 
(wetting)  phase  of  water  uptake.  During  the  drying  phase,  Do  is  found  to  be 
typically  larger  by  a  factor  of  10,  and  logio(Do)  executes  a  sharp  increase  for  (dry¬ 
ing)  9  values  between  50  >6  (kg/m3)  >  0.0.  These  features  are  consistent  with 
the  aforementioned  hysteresis  of  0,  and  have  been  attributed  by  Pel  (1995)  to  the 
transition  from  liquid  to  vapor- dominated  moisture  transport. 


Moisture  content  (grams/cnv  ) 


Figure  2.  Diffusion  coefficient  variation  with  moisture  content,  for  several  common  porous 
building  materials  (van  Zijl  1999). 
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Thus  while  De  performs  a  pivotal  function  in  Equation  1,  its  ultimate  functional¬ 
ity  is  constrained  by  the  isothermal  sorption  curve  (Figure  1).  This  is  particu¬ 
larly  relevant  when  the  RH  is  changing.  As  the  water  vapor  diffuses  into  the 
material  according  to  the  diffusion  coefficient,  it  is  being  absorbed  by  the  pore 
wall  structure,  according  to  the  sorption  coefficient.  A  material  with  both  a  high 
diffusivity  and  a  high  absorption  coefficient  will  appear  to  retard  the  transport  of 
water  vapor  under  dynamic  circumstances  because  the  water  is  absorbed  before 
it  can  appreciably  diffuse. 

In  recent  years,  percolation  theory  has  offered  a  valuable  paradigm  for  gleaning 
insights  into  the  defining  physics  responsible  for  many  of  these  features;  the  re¬ 
view  by  Celia  et  al.  (1995),  is  particularly  helpful.  In  addition,  the  advent  of  new 
technology,  such  as  NMR,  has  enabled  the  detailed  measurements  of  actual  in 
situ  moisture  profiles  in  real  time  experiments.  These  observations  provide  the 
crucial  data  for  the  sorption  and  desorption  of  moisture  coefficients  under  a  vari¬ 
ety  of  environmental  conditions  and  materials.  The  reader  is  referred  to  the  ex¬ 
cellent  works  of  Pel  (undated). 

Finally  for  realistic  modeling  of  any  physical  system,  one  must  be  aware  of  the 
impact  that  any  ancillary  boundary  conditions  can  have  on  a  set  of  (perceived) 
primary  boundary  conditions.  This  constraint  can  be  particularly  crucial  for  a 
porous  system  containing  both  invading  and  receding  fluids  (Aker  1996;  Des¬ 
camps  1997).  Descamps  has  observed  that  water  imbibition  into  porous  building 
materials  is  a  two-phase  flow  process  in  which  the  imbibing  water  displaces  the 
air  that  is  initially  present  in  the  pore  space.  Accurate  non-destructive  meas¬ 
urements  of  transient  moisture  profiles  during  imbibition  into  calcium  silicate 
brick  indicate  that  the  air  outflow  boundary  condition  has  a  predominant  influ¬ 
ence  on  the  water  imbibition  process.  As  a  consequence,  significant  differences 
have  been  observed  in  the  water  sorption  and  transport  coefficients  derived  from 
water  imbibition  experiments  with  free  and  bounded  air  outflow  boundary  condi¬ 
tions. 


Diffusion  of  Moisture  in  Polymers 

The  set  of  Equations  (1  to  3)  and  the  attendant  diffusivity  functions,  defines  the 
basic  diffusion  mechanisms  for  moisture  migration  through  the  masonry  portion 
of  the  masonry-FRP  applique.  The  same  equation  set  exists  for  each  region,  but 
with  a  different  De  and  Dt  functionality  and  attendant  sorption  curves.  In  par¬ 
ticular,  the  literature  concerning  the  effects  of  fluids  on  polymeric  composites 
has  increased  dramatically  in  recent  years,  in  part  due  to  the  immense  variety  of 
materials  and  circumstances  associated  with  these  compounds.  Polymers  are 
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highly  complex  materials  that  vary  in  structure  and  physico-chemical  properties. 
FRP  composites  adjoin  an  assortment  of  fibers  with  intricate  fiber/matrix  inter¬ 
phases  and  bondings,  further  complicating  the  resultant  composite.  The  com¬ 
plexity  is  further  confounded  by  the  sorption  of  many  different  fluids,  each  with 
its  own  unique  interaction  between  the  polymer,  fiber,  and  interface  within  the 
composite  (Weitsman  1995). 

The  primary  issues  associated  with  the  mechanical  performance  of  FRP  compos¬ 
ites  in  the  presence  of  fluids,  concerns  their  dimensional  stability,  strength,  fa¬ 
tigue  response,  and  impact  resistance  above  and  below  the  glass  transition  tem¬ 
perature,  Tg.  To  a  large  degree,  these  issues  are  still  only  partially  understood, 
even  in  the  absence  of  fluids.  Polymers  are  somewhat  porous,  but  share  many  of 
the  diffusion  and  sorption  considerations  of  their  ceramic  counterparts,  including 
hysteresis.  However  unlike  the  typical  (cured)  ceramic,  polymeric  volume 
changes  can  be  quite  sensitive  to  the  infusion  of  a  variety  of  penetrant  liquids. 
Typical  volume  strains  of  2  percent  are  common  for  comparable  water  weight- 
percent  uptake  in  a  sorptive  polymer.  Finally,  the  application  geometry  of  the 
polymer  can  play  a  key  role  in  the  effective  diffusion  coefficient,  and  increases  of 
several  orders  of  magnitude  have  been  reported  for  polymeric  films  in  going  from 
a  dry  to  moist  state  (Valen  1998). 

To  quantify  the  many  features  found  in  the  typical  polymer  system  undergoing 
fluid  uptake  (i.e.,  moisture),  Lefebvre  et  al.  (1989)  set  forth  a  model  for  the  non¬ 
linear  diffusion  coefficient  formulated  on  a  methodology  to  relate  the  diffusion 
coefficient  of  small  penetrant  molecules  to  temperature,  strain,  and  penetrant 
concentration.  The  approach  is  based  on  well  known  free  volume  theories  of 
Doolittle  and  others  (Ferry  1980),  and  assumes  that  the  transport  kinetics  are 
governed  by  the  constant  redistribution  of  the  free  volume,  caused  by  the  seg¬ 
mental  motions  of  the  polymeric  chains.  In  essence,  this  strain  induced  volume 
fluctuation  mimics  a  pumping  action  that  augments  the  normal  diffusion  proc¬ 
ess,  which  requires  the  penetrant  molecule  to  surmount  an  energy  barrier  re¬ 
lated  to  the  fluctuating  pore  size.  Their  expression  for  the  diffusion  coefficient, 
De  below  the  glass  transition  temperature  Tg,  is: 

De  =  D0exp(-AH/RT)exp(A/B),  Eq  5 

in  which: 


A  =  (Polf0)(t  +  yGt4),  and  B  =  (f0  +  £+y6N) 
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where: 

fo  =  the  fractional  free  volume  of  the  polymer,  equal  to  the  void  volume  di¬ 

vided  by  the  specific  volume 

s  =  volume  dilatation  of  the  free  volume  due  to  external  loads 

G  =  the  penetrant  concentration,  in  this  case  moisture 

Y  -  the  volumetric  coefficient  of  swelling,  as  the  result  of  the  penetrant  ad¬ 

sorption 

N  =  an  adjustable  parameter  related  to  the  non-zero  volume  change  occurring 
during  moisture  uptake 

Bd  =  a  numerical  parameter  inversely  related  to  the  critical  pore  size,  and  R  is 
the  gas  constant 

AH  =  an  activation  energy  related  to  the  effective  potential  barrier  that  must 
be  surmounted  by  the  penetrant  molecule  to  percolate  through  the 
critical  pore  size 

T  =  the  absolute  thermodynamic  temperature,  measured  in  degrees  Kelvin 

Do  —  the  diffusion  coefficient  in  the  reference  state. 


Equation  5  indicates  that  the  free  volume  continues  to  govern  the  effect  of  stress 
and  swelling  below  the  glass  transition  temperature,  Tg,  where  the  fractional 
volume  is  much  less  temperature  dependent.  Below  Tg,  this  metastable  (super¬ 
cooled)  state  causes  the  free  volume  to  collapse  slowly  with  time  until  an  equilib¬ 
rium  is  reached.  The  phenomena  is  known  as  physical  aging.  Any  successful 
diffusion  modeling  in  the  glassy  range  must  incorporate  this  important  effect. 
The  implication  is  that  these  forces  act  primarily  on  very  local  volumetric  prop¬ 
erties  such  as  the  critical  void  size. 


By  contrast,  temperature  appears  to  affect  motions  over  larger  volumes  encom¬ 
passing  cooperative  segmental  motions.  The  Arrhenius  temperature  dependence 
shown  in  Equation  5  has  been  observed  experimentally  (Comyn  1985)  in  numer¬ 
ous  penetrant-polymer  systems  below  Tg.  Further,  the  Arrhenius  form  may  also 
be  used  above  Tg,  with  an  activation  energy  (AH)  differing  from  that  below  Tg, 
which  accounts  for  the  abrupt  change  in  the  temperature  dependence  at  the 
glass  transition.  Lefebvre  et  al.  (1989)  proceed  by  explaining  that,  for  tempera¬ 
tures  above  Tg,  the  e  terms  in  A  and  B  of  Equation  5  acquire  an  additional  tem¬ 
perature  sensitivity  proportional  to  ( T -  7^).  The  applicability  of  the  Arrhenius 
form  is  based  in  thermodynamic  phase  transitions  and  is  further  justified  on  the 
basis  of  the  analytical  suppleness  of  the  actual  physical  data. 

In  polymeric  materials,  Equation  5  affords  an  efficient  parameterization  of  De  in 
terms  of  parameters  with  a  readily  identifiable  physical  interpretation.  As  in 
the  case  of  ceramics,  the  most  reliable  values  for  these  parameters  are  rendered 
from  actual  experimental  data.  Table  1  (Lefebvre  et  al.  1989)  lists  some  repre¬ 
sentative  values  for  T  >  Tg. 
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Table  1.  Some  experimental  values  of  fitting  parameters  used  to  define  Dg  in 
Equation  5,  for  several  listed  polymeric  materials. 


Fitting-  Parameter 

Ethylbenzene- 

Polystyrene 

Toluene- 

Polystyrene 

Toluene- 

Polyvinylacetate 

fiolfo 

17.9 

16.4 

8.8 

fo 

0.074 

0.014 

0.319 

Ta(C<) 

80 

88 

50.3 

N 

0.61 

0.90 

0.91 

r 

0.48 

0.51 

2.10 

Logi0(Do[m2/s]) 

-24.3 

-24.6 

-15.5 

However,  to  demonstrate  the  predictive  capability  of  their  model,  this  work  com¬ 
pared  these  experimental  fitting  parameters  with  the  expected  constants  deriv¬ 
able  from  the  model’s  rudimentary  beginnings  in  Wells-Landel-Ferry  theory 
(Ferry  1980).  The  values  obtained  for  Brifo,  Do  and  Tg  were  in  reasonable  accord 
with  empirical  data,  and  found  to  be  independent  of  the  nature  of  the  solvent. 
Some  values  for  the  enthalpic  activation  energy,  A H,  derived  from  the  literature 
are:  21.8  kJ/mol  for  Vinyl  Ester  water  diffusion,  69.74  kJ/mol  for  Epoxy  resin 
water  diffusion,  and  33.23  kJ/mol  for  the  diffusion  of  water  in  Isopolyester  (Chin 
et  al.  1999). 
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3  Computer  Modeling  and  Simulation 


The  typical  unreinforced  masonry  wall  requiring  an  FRP  upgrade  is  located  on 
the  exterior  surface  of  a  building,  and  is  a  load-bearing  structure  of  the  building 
envelope.  Masonry  walls  located  within  the  interior  may  also  require  some  rein¬ 
forcement  against  seismic  events,  but  typically  are  not  load  bearing.  This  study 
assumed  that  the  exterior  wall  is  constructed  of  a  single  course  of  CMUs  without 
external  brick  facing.  As  such,  the  large  majority  of  FRP  appliques  will  be  ap¬ 
plied  to  the  interior  surface,  where  the  environmental  conditions  of  moisture  and 
humidity  are  more  stable  and  fluctuations  less  extreme.  If  the  interior  is  kept 
within  a  comfortable  range  for  (say)  office  environments  (30  to  60  percent  RH, 
and  72  to  78  °F),  then  common  architectural  and  construction  guidance  indicates 
that  moisture  will  flow  towards  the  cooler  and  dryer  surface.  Furthermore,  labo¬ 
ratory  studies  have  shown  that,  for  the  sake  of  both  economy  and  mechanical  re¬ 
sponse,  FRP  reinforcement  in  the  form  of  fabric  strips  is  preferable  to  fabric  that 
covers  the  whole  surface  of  masonry  walls  (Triantafillou  1998;  Rivera  and 
Karbhari  2003). 

Drawing  on  these  facts,  an  initial  simulation  modeling  effort  was  made  using  a 
wall  built  from  CMU.  Figure  3  shows  a  test  section  of  such  a  wall  along  with  the 
dimensions  (in  mm)  of  the  individual  CMUs.  The  grid  structure  that  overlays 
the  individual  CMU  (viewed  here  from  the  vertical)  is  indicative  of  the  mesh  de¬ 
velopment  required  for  a  finite  element  calculation.  Note  that,  in  this  figure,  the 
FRP-applique  appears  to  cover  the  entire  interior  surface  of  the  CMU  when  in¬ 
stalled  in  the  traditional  vertical  position.  Since  total  coverage  of  the  wall’s  sur¬ 
face  is  not  recommended,  for  simulation  a  coverage  of  about  60  percent — an  ex¬ 
treme  upper  limit  of  FRP  usage — will  be  used. 

Figure  3b  shows  an  important  feature  of  a  typical  CMU.  The  CMU  has  large 
empty  bays  (voids)  that  can  readily  communicate  with  the  surrounding  internal 
thermal  and  moisture  environment  from  above  and  below.  In  actual  usage, 
these  internal  voids  can  (and  do)  have  long  indirect  paths  that  have  contact  with 
some  unknown  combination  of  interior  and  exterior  air  conditions.  Given  the 
requirement  of  well  defined  boundary  conditions  for  the  solution  of  the  combined 
Equation  set  1  to  3  above,  this  work  considers  a  very  restrictive  set  of  thermal, 
moisture,  and  material  thickness. 
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b. 

Figure  3.  Test  section  of  a  wall  (Figure  3a)  constructed  from  CMUs,  along  with  an  individual  CMU’s 
dimensions.  The  2  mm  of  FRP  composite  thickness  (Figure  3b),  has  been  exaggerated  for  clarity  at  this 
scale. 
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Thus,  for  simulation  purposes,  the  ~36  mm  of  solid  CMU  portion  nearest  the 
FRP  applique,  will  be  considered  as  the  defining  relevant  thickness  of  concrete, 
and  a  corresponding  FRP  thickness  of  ~2.1  mm.  Results  derived  from  these 
simulations  must  be  viewed  as  somewhat  extreme  so  that  actual  conditions  of 
moisture  and  temperature  gradients  will  be  less  severe. 

The  final  computational  simulation  is  then  reduced  to  that  of  modeling  a  single 
layer  of  solid  CMU  in  natural  contact  with  an  FRP  applique,  that: 

3.  Is  approximately  0.1  to  0.066  the  thickness  of  the  CMU  base  (Figure  4a). 

4.  Has  the  initial  interior  surface  constrained  to  a  hygrothermal  environmental 
range  of  30  to  60  percent  RH,  and  72  to  78  °F,  with  an  exterior  surface  where  the 
ambient  outside  humidity  and  temperature  are  free  parameters,  houtl  Tout,  re¬ 
spectively  (Figure  5a). 

5.  Has  a  60  percent  lateral  surface  coverage,  of  FRP  applique. 

6.  Is  surrounded  by  repeating  and  identical  CMU-FRP  units  in  the  “x”  and  “z”  di¬ 
rection  of  Figure  4a.  This  is  also  known  as  a  “periodic  boundary  condition,” 
therefore  the  only  thermal  and  moisture  flux  is  in  the  “X-Y”  direction.  An  actual 
three-dimensional  model  is  needed  for  combined  hygrothermal-mechanical 
analysis  (Figure  4b).  However  under  the  consideration  of  these  boundary  condi¬ 
tions,  the  mass  and  heat  transport  problem  can  be  reduced  to  a  two-dimensional 
model,  subject  to  combined  temperature  and  moisture  gradients. 

Apart  from  the  symmetry,  this  is  required  to  reduce  the  scale  of  numerical  computa¬ 
tion  and  the  anticipated  time-consuming  nonlinear  numerical  procedures.  See  Fig¬ 
ure  5a,  for  details  and  the  subsequent  reduction  of  the  3D  model  to  the  corresponding 
2D  analog.  No  attempt  is  made  to  account  for  or  to  differentiate  the  mortar  joints 
connecting  the  individual  CMU-FRP  units  found  in  a  traditional  layered  wall. 


Numerical  Solutions  and  Finite  Element  Implementation 

Given  the  geometry  and  the  boundary  conditions  required  for  modeling  the 
structures  in  Figure  4,  one  can  now  attempt  to  unravel  Equation  set  1  to  3  for  a 
solution.  General  analytic  solutions  involving  this  set  of  coupled,  nonlinear  dif¬ 
ferential  equations  are  not  possible,  even  for  restricted  ID  and  2D  geometry. 
The  problem  is  further  confounded  by  the  use  of  actual  physical  data  (Appendix 
A)  that  is  functionally  dependent  on  the  field  variables  of  moisture,  6,  and  tem¬ 
perature,  T,  and  any  subsequent  structural  analysis  are  by  their  nature 
three-dimensional  problems. 
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Figure  4.  3D  model  of  the  FRP  reinforced  concrete  block:  (a)  schematic  diagram  and  (b)  FE 
model  mesh  (Dimensions  are  unitless,  meant  to  indicate  proportion  only). 
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b. 

Figure  5.  Reduced  2D  model  of  the  original  3D  problem  used  for  analysis  of  thermal  and 
moisture  transfer.  Figure  (5a)  depicts  lateral,  interior  and  exterior  boundary  conditions;  (5b) 
represents  corresponding  FE  model  mesh. 


20 


ERDC/CERL  TR-03-20 


Fortunately,  this  set  of  equations  is  responsive  to  general  numerical  analyses 
and  solution,  in  particular  that  of  the  finite  element  method.  The  finite  element 
numerical  method  provides  great  flexibility  for  complex  model  geometry  and  al¬ 
lows  for  parametric  studies  of  various  FRP  reinforcement  arrangements,  which 
would  be  very  difficult  if  not  impossible  using  traditional  theoretical  analysis.  In 
addition,  the  actual  level  of  numerical  modeling  effort  required  for  the  detailed 
solution  of  these  equations,  and  their  attendant  boundary  conditions,  is  non¬ 
trivial.  Numerical  modeling  services  for  this  report  were  obtained  through  a 
contractual  agreement  between  ERDC-CERL  and  Clark  Atlanta  University.^ 


A  number  of  finite  element  packages  are  commercially  available.  This  work  used 
the  ANSYS  commercial  package  in  some  of  the  early  preliminary  analysis.  How¬ 
ever,  the  ANSYS  package  does  not  have  the  capability  of  performing  coupled  hy- 
grothermal-mechanical  analysis.  (In  fact,  no  commercial  finite  element  software 
packages  have  this  capability.)  To  implement  the  finite  element  analysis  proce¬ 
dure,  the  FEMLAB  finite  element  analysis  software  was  chosen  for  the  develop¬ 
ment  of  the  coupled  hygrothermal-mechanical  analysis  procedure.  The  FEMLAB 
software  is  built  on  the  MATLAB  mathematical  tool  using  matrix  operation.  The 
program  is  primarily  designed  for  solving  differential  equations  using  the  finite 
element  method.  The  FEMLAB  finite  element  package  allows  users  to  modify 
the  predefined  finite  element  formulation  and  to  input  any  type  of  partial  differ¬ 
ential  equations  to  be  solved  in  a  defined  geometry.  The  coefficients  of  the  equa¬ 
tions  generally  can  be  functions  of  independent  and  dependent  variables. 


Replacing  the  moisture  variable  with  its  equivalent  representation  in  terms  of 
humidity,  h,  yields  the  Galerkin  finite  element  formulation  of  the  coupled  hu¬ 
midity  and  heat  transport  equations: 


Eq  6 


where: 

V  =  the  volume  of  a  finite  element 

S  =  surface  where  the  vector  of  humidity/heat  flux  is  specified 

q*  =  boundary  humidity/heat  flux 


^  Lin,  Mark  W.,  Department  of  Mechanical  and  Aerospace  Engineering,  University  of  Alabama  in  Huntsville,  AL 
35899;  The  prime  contractor  was  Department  of  Engineering  and  Chemistry,  Clark  Atlanta  University,  Atlanta,  GA 
30314.  All  inquiries  concerning  the  numerical  details  of  this  report  should  be  addressed  to  Dr.  Lin,  at  the  above 
address;  tel:  (256)  824  6325;  or  e  mail:  lin@mae.uah.edu. 
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the  vector  of  an  allowable  virtual  humidity/temperature 


W= 


dx  dy  dz 

0  0  0 


0  0 


0 

_d_  d_  d_ 
dx  dy  dz  J 


is  a  differential  operator. 


In  the  above  equation, 


C 

\D* 


is  the  element  humidity/heat  capacity  matrix 
is  the  element  humidity/heat  diffusion  matrix. 


Based  on  the  coupled  humidity/temperature  governing  equations,  Equations  1,  2, 
and  3: 

[c*  ]  and  [d*  ]  have  the  following  matrix  forms: 


0 

pc 


and 


0 

0 

0 

0 

0 


0  0  DTx  0  0 

Dhy  0  0  DTy  0 

0  Dhz  0  0  DTz 

OO^OO 
OOO^O 
0  0  0  0  4 


Using  element  shape  functions  defined  in  the  matrix  form  as 


{V(x,y,z)} 


JX  0  n2  0 
[0  X  0  n2 


where  r  is  the  total  number  of  nodes  per  element,  and  the  element  humid¬ 
ity/temperature  vector  can  be  discretized  in  terms  of  nodal  humid¬ 
ity/temperature  vector: 
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(i  =  l,2,...,r) 


as: 


Eq  7 


Substituting  Equation  7  into  Equation  6,  the  finite  element  equation  becomes: 

dh/ 

r  h, 


[CJ 


dt 

ST 


l  dt) 


+[D'J L  ={&} 


Eq  8 


where: 


[Dj]=f  [B]t[D'\BW  —  element  humidity/heat  diffusion  matrix 

JV 

[c:]=\\nIc-\nUv  —  element  humidity/heat  capacity  matrix 
{Qe}=  j  IN}#  dS—  element  boundary  humidity/heat  flux 

[i]=W 


In  the  above  finite  element  equation,  the  shape  functions,  {N(x,y,z)\,  are  poly¬ 
nomials  in  terms  of  the  local  coordinates  of  a  finite  element  and  are  used  to  ap¬ 
proximate  the  variations  of  the  state  variables,  i.e.,  humidity  and  temperature, 
within  each  element.  The  values  of  the  state  variables  at  any  location  within  the 
element  can  then  be  evaluated  based  on  the  nodal  values  of  the  state  variables 
as  expressed  in  Equation  7.  The  order  of  the  polynomial  used  in  the  shape  func¬ 
tion  will  dictate  the  accuracy  and  the  characteristics  of  the  finite  element  nu¬ 
merical  solution.  In  general,  the  higher  the  order  of  the  shape  function,  the 
higher  the  accuracy  of  the  solution.  Nevertheless,  it  requires  larger  number  of 
nodes  per  element  and  longer  computation  time  to  reach  a  solution. 


Higher  order  nonlinear  shape  functions  are  generally  more  suitable  for  models 
with  curved  boundaries.  Since  there  are  no  curved  boundaries  in  the  model  em¬ 
ployed,  linear  shape  functions  were  used  in  the  present  analysis.  Such  a  choice 
was  based  on  the  considerations  for  reducing  computation  time  of  the  analysis  in 
three-dimensional  models  as  well  as  for  avoiding  the  need  to  account  for  nonlin¬ 
earities  both  in  material  properties  and  model  shape  functions. 

To  solve  the  transient  problem,  the  finite  element  equation,  Equation  8,  was  dis¬ 
cretized  by  the  generalized  trapezoidal  rule  with  respect  to  the  time  coordinate. 
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The  resulting  system  of  equations  at  the  time  instant  was  then  solved  by  the  it¬ 
erative  linearization  method.  Note  that,  at  the  first  time  instant,  the  diffusion 
matrix  and  the  capacity  matrix  were  computed  based  on  the  initial  humidity  and 
temperature  conditions.  Subsequently,  their  values  were  updated  based  on  the 
solution  of  the  humidity  and  temperature  fields  obtained.  These  values  were 
then  used  for  the  next  time  step.  The  new  system  of  equations  was  solved  again 
and  the  procedure  continued  with  increasing  time  instant  until  it  reached  the 
final  desired  time  instant. 


Numerical  Solutions  for  Moisture,  Total  FRP-CMU  Coverage 

With  the  above  finite  element  modeling  tools  at  hand,  a  simple  test  case  was  run 
to  provide  a  cross  check  for  the  proper  functioning  of  the  coupled  computer  codes, 
and  a  reality  check  from  both  physical  and  field  experience.  For  this  elementary 
case,  the  2D  problem  shown  in  Figure  5  was  further  reduced  to  a  ID  problem  by 
having  the  FRP  totally  cover  the  CMU.  The  interior  FRP  surface  was  con¬ 
strained  to  50  percent  RH  and  20  °C,  while  the  exterior  surface  RH  was  held  at 
97  percent  RH  with  the  temperature  permitted  to  vary  from  10  to  60  °C.  As 
such,  this  simulation  was  representative  of  extreme  hygrothermal  external  envi¬ 
ronments.  Figures  6  and  7,  respectively,  show  spatial  distribution  solutions  (in  a 
direction  normal  to  FRP-CMU  interface)  of  humidity,  h,  and  total  moisture,  0. 
Each  Figure  is  a  composite  of  both  FRP  and  CMU  steady-state  data  for  several 
temperatures  extending  from  10  to  60  °C.  For  this  project,  the  FRP  has  a  stan¬ 
dard  thickness  of  1.94  mm  and  the  CMU  base  has  a  standard  thickness  of  19.4 
mm. 

Figure  6  shows  the  humidity  starting  at  50  percent  and  rapidly  rising  monotoni- 
cally  towards  the  FRP-CMU  interface,  where  there  is  a  spread  in  RH,  h,  associ¬ 
ated  with  temperature  variations.  The  relative  humidity  fulfills  the  humidity 
continuity  boundary  condition  (Equation  2)  at  the  interface,  and  has  a  relatively 
smaller  spread  as  it  tends  toward  the  97  percent  RH  boundary  condition  on  the 
outer  surface  of  the  CMU.  This  figure  confirms  the  physical  notion  that  most  of 
the  RH  drop  through  the  system  will  occur  in  the  (oil-based)  FRP.  The  CMU  is 
relatively  porous  and,  by  comparison,  more  transparent  to  moisture  migration. 
Thus  the  RH  will  be  more  or  less  constant  across  the  CMU. 

Further  inspection  reveals  that  the  RH  at  the  FRP-CMU  interface  increases 
monotonically  with  the  outside  temperature.  Mathematically,  this  arises  from 
the  DfJT,  term  in  Equation  1,  and  again  confirms  the  physical  experience  that  a 
temperature  gradient  can  act  as  a  pseudo-moisture  source. 
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X  (Concrete,  units  of  0.1  mm) 
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X  (FRP,  units  of  0.1  mm) 


Figure  6.  ID  steady-state  moisture  transport  problem,  humidity  profile,  with  FRP  totally 
covering  the  CMU;  the  interior  FRP  surface  is  constrained  to  50%  RH  and  20  °C,  while  the 
exterior  surface  RH  is  held  at  97%  RH  with  the  temperature  permitted  to  vary  from  10  to  60  °C. 


Figure  7.  1 D  steady-state  moisture  transport  problem,  total  moisture  density  0  profile,  with  FRP 
totally  covering  the  CMU;  the  interior  FRP  surface  is  constrained  to  50%  RH  and  20  °C,  while  the 
exterior  surface  RH  is  held  at  97%  RH  with  the  temperature  permitted  to  vary  from  10  to  60  °C. 
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Closer  inspection  however,  reveals  a  cautionary  feature.  For  temperatures 
greater  than  45  °C,  it  is  seen  that  the  RH  is  greater  than  100  percent.  Such  su¬ 
persaturation  might  suggest  pooling  or  condensation  of  moisture  in  a  laboratory 
experiment.  However  there  is  no  physical  mechanism  for  this  condensation  in¬ 
corporated  into  the  isothermal  absorption  curves  for  either  FRP  or  CMU  (Figure 
1).  This  anomaly  appears  to  originate  from  the  fact  that  the  isothermal  absorp¬ 
tion  curves  for  FRP  and  CMU  (Appendix  A)  are  analytically  continuous  for  all 
values  of  the  relative  humidity,  h,  even  though  the  physical  region  of  applicabil¬ 
ity  is  constrained  to  a  RH  of  98  percent  or  less  (see  the  end  of  section  2.1).  Thus 
the  computer  is  able  to  converge  a  solution  of  questionable  physical  reality. 
However  the  trend  is  one  that  has  a  basis  in  reality.  The  thermal  gradient  is  al¬ 
ways  from  a  warm  surface  toward  a  colder  surface.  As  such,  there  exists  a  clear 
tendency  for  increased  moisture  at  the  cooler  interface.  An  immediate  conse¬ 
quence  is  to  avoid  total  coverage  of  CMU  in  such  extremes  of  humidity  and  tem¬ 
perature.  This  conclusion  derives  further  support  from  the  experience  of  exterior 
paint  applications  in  cold  climes  (Korhonen  and  Bayer  1989),  which  cautions  to 
provide  ample  ventilation  of  painted  exterior  masonry,  to  avoid  paint  delamina¬ 
tions,  primarily  from  efflorescence. 

Figure  7,  by  comparison,  has  many  of  the  qualitative  features  found  in  Figure  6 
with  one  significant  difference,  the  discontinuity  of  the  moisture  density,  6, 
across  the  FRP-CMU  interface.  This  anticipated  difference  in  moisture  content 
makes  (petroleum-based)  paints  effective  in  preserving  wood  from  moisture  up¬ 
take.  Furthermore,  a  detailed  analyses  of  the  interfacial  moisture  gradients,  the 
q  +  D$J6  +  D'NT-  0  of  Equation  1,  certifies  that  the  moisture  flux  requirement 
of  Equation  2  is  also  fulfilled.  In  summary,  these  collective  results  support  the 
conclusion  that  the  finite  element  computer  program  can  provide  useful  results 
in  keeping  with  physical  intuition  and  field  data,  provided  that  the  calculation 
returns  values  of  the  relative  humidity  within  the  hygroscopic  range  of  a  given 
isothermal  sorption  curve. 


Numerical  Solutions  for  Moisture,  Partial  FRP-CMU  Coverage 

Drawing  on  the  success  of  the  previous  ID  model,  a  calculation  was  performed 
for  the  2D  case  shown  in  Figure  5.  In  this  rendering,  the  FRP  covers  about  60 
percent  of  the  interior  lateral  surface  (X-Z,  plane).  As  in  the  ID  case,  the  interior 
FRP  surface  was  constrained  to  50  percent  RH  and  20  °C,  while  the  exterior  sur¬ 
face  RH  was  held  at  97  percent  RH  with  the  temperature  permitted  to  vary  from 
10  to  60  °C.  Spatial  distribution  solutions  of  humidity,  h,  and  total  moisture,  Q, 
were  computed  throughout  the  2D  structure  (within  the  Y-X  plane).  However 
the  middle  of  the  CMU-FRP  block  contains  a  plane  of  symmetry  with  respect  to 
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the  “X”  direction.  Calculations  performed  in  this  plane  (containing  a  “Y”  axes) 
will  have  the  highest  values  of  humidity,  h,  and  total  moisture,  9,  relative  to  the 
corresponding  uncovered  portions  of  CMU.  As  such  these  centerline  values  pro¬ 
vide  a  reasonable  benchmark  for  comparison  between  the  ID  and  2D  cases. 

Figures  8  and  9,  respectively,  show  the  steady-state  humidity  and  total  moisture 
as  function  of  spatial  distribution.  Each  Figure  is  a  composite  of  both  FRP  and 
CMU  data  for  several  temperatures  extending  from  10  to  60  °C,  and  have  the 
same  FRP-base  dimensions  as  before.  Many  of  the  global  physical  characteris¬ 
tics  of  these  Figures  are  analogous  to  the  ID  case.  However,  the  relative  humid¬ 
ity  remains  within  physical  limits  for  all  temperatures.  As  in  the  ID  case,  there 
is  an  analogous  increase  of  humidity,  h,  and  total  moisture,  9,  with  increasing 
external  temperature  arising  from  the  D'NT term  in  Equation  1.  Further  inspec¬ 
tion  of  Figures  8  and  9  reveals  that  both  the  humidity  and  total  moisture  are  ac¬ 
tually  less  at  the  FRP-CMU  interface,  than  at  the  external  boundary.  This  is  the 
case  for  external  temperatures  of  50  °C  or  less,  and  again  underscores  the  virtue 
of  ventilating  the  CMU,  even  with  only  a  minor  fraction  of  the  total  surface  area. 

While  Figures  8  and  9  provide  quantitative  information  within  the  symmetry 
plane,  Figure  10  shows  a  more  dramatic  portrayal  of  the  overall  moisture  flux 
(within  the  X-Y  plane).  Figure  10  clearly  depicts  the  evacuation  routes  available 
to  moisture  that  might  otherwise  be  trapped  behind  the  CMU-FRP  interface. 
Figure  1 1  shows  the  resulting  combination  if  the  X-Y  humidity  profile  is  super¬ 
imposed  on  Figure  10.  The  moisture  flux  shown  in  Figure  11  appears  to  be  nor¬ 
mal  to  the  surfaces  of  constant  humidity,  which  is  an  expected  result  given  that 
the  flux,  q,  scales  proportional  to  V0. 

Finally,  additional  information  consolidated  within  Figures  8  and  9  relates  to  the 
thermal  profiles.  To  converge  a  consistent  solution  for  Equation  set  1  to  3  involv¬ 
ing  the  total  moisture,  9,  an  equivalent  set  of  (steady-state)  thermal  information 
is  also  generated,  T-  T\r,  9).  Restricting  this  solution  set  to  the  symmetry  plane 
implied  for  Figures  8  and  9  results  in  the  thermal  profile  of  Figure  12.  As  in 
those  cases,  the  Figure  is  a  composite  of  both  FRP  and  CMU  steady-state  data, 
with  the  interior  FRP  surface  constrained  to  50  percent  RH  and  20  °C  (Centi¬ 
grade),  while  the  exterior  surface  RH  was  held  at  97  percent  RH  with  the  tem¬ 
perature  permitted  to  vary  from  10  to  60  °C. 
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X  (  Concrete,  units  of  0.1  mm) 


0  38.8  77.6  116.4  155.2  194 


X  (FRP,  units  of  0.1  mm) 


Figure  8.  2D  steady-state  moisture  transport  problem,  humidity  profile,  with  FRP 
partially  covering  (-60%)  the  CMU;  data  taken  from  the  (Y-Z)  symmetry  plane;  the  interior 
FRP  surface  is  constrained  to  50%  RH  and  20  °C;  while  exterior  surface  RH  is  held  at 
97%  RH  with  the  temperature  permitted  to  vary  from  10  to  60  °C. 


Figure  9.  2D  steady-state  moisture  transport  problem,  total  moisture  density  0  profile, 
with  FRP  partially  (-60%)  covering  the  CMU;  data  taken  from  the  (Y-Z)  symmetry  plane; 
the  interior  FRP  surface  is  constrained  to  50%  RH  and  20  °C,  while  the  exterior  surface 
RH  is  held  at  97%  RH  with  the  temperature  permitted  to  vary  from  10  to  60  °C. 
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Figure  10.  2D  steady-state  moisture  transport  problem,  depicting  total  moisture 
flux  profile  (within  the  X-Y  plane),  and  FRP  partially  (~60%)  covering  the  CMU;  this 
chart  delineates  evacuation  routes  available  to  moisture  that  might  otherwise  be 
trapped  behind  the  CMU-FRP  interface;  the  interior  FRP  surface  is  constrained  to 
50%  RH  and  20  °C;  the  exterior  surface  RH  is  held  at  97%  RH  with  a  temperature 
of  order  20  °C. 


Figure  11.  Figure  10,  with  the  corresponding  X-Y  humidity  profile  superimposed. 
Note  that  the  moisture  flux  appears  to  be  normal  to  the  surfaces  of  constant 
humidity,  which  is  an  expected  result  given  that  the  flux,  q,  scales  proportional  to 
V0,  (see  Equation  1). 
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X  (Concrete,  units  of  0.1  mm) 

0  38.8  77.6  116.4  155.2  194 


— B— T=60  C 
T=50  C 
— *-T=40  C 
-b-T=30  C 
“ T=20  C 
— T=10  C 


0  3.88  7.76  11.64  15.52  19.4 

X  (FRP,  units  0.1  mm) 


Figure  12.  FRP-CMU  thermal  profile,  with  FRP  partially  (-60%)  covering  the  CMU; 
data  taken  from  the  (Y-Z)  symmetry  plane;  the  interior  FRP  surface  constrained  to 
50%  RH  and  20  °C;  exterior  surface  RH  is  held  at  97%  RH  with  the  temperature 
permitted  to  vary  from  10  to  60  °C. 

Figure  12  shows  an  almost  linear  profile  of  temperature  as  a  function  of  distance 
through  each  of  the  FRP-CMU  structures.  This  quasi-linear  profile  is  the  signa¬ 
ture  of  an  essentially  constant  thermal  conductivity  parameters,  kr,  for  the  FRP- 
CMU  system.  A  simple  thermal  flux  calculation  based  on  constant  FRP-CMU 
material  parameters  (Appendix  A),  yields  an  interfacial  temperature  quite  close 
to  that  resulting  from  the  more  exacting  calculation  shown  in  Figure  12.  This 
observation  verifies  the  prediction  made  at  the  end  of  section  2.3,  regarding  the 
role  of  temperature  as  a  higher  order  parameter  for  the  convergence  of  the  total 
moisture  solution.  The  resulting  insensitivity  tends  to  decouple  the  two  solu¬ 
tions,  and  makes  for  efficient  computer  coding,  at  least  for  these  particular  build¬ 
ing  materials. 

Closer  inspection  reveals  a  slight  curvature  in  the  thermal  profiles  of  concrete 
and  FRP.  This  reflects  the  small  variation  in,  kr,  due  to  any  explicit  moisture 
sensitivity,  as  well  the  very  slight  contributions  arising  from  the  V(hPs(  T)  terms 
found  in  the  “WUFI-ORNL/IBF’  formulation  (Kuenzel  et  al.  undated  [Oak 
Ridge])  of  DFJT( see  Appendix  A,  Equation  A5).  In  fact,  a  detailed  analyses  of 
the  thermal  flux  vector: 


FRP  Concrete 

30 

20 


qTh  =  {kT  +  hv(6p(7)///)(AH//?7£)(/iPs(7))}vr 
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using  the  steady-state  numerical  data  of  Figures  12  and  9  reveals  that  the  vector 
is  conserved  throughout  the  symmetry  plane  of  the  CMU-FRP.  A  consequence  of 
this  insensitivity  to  latent  heat  terms  is  that  the  interfacial  temperature  with 
and  without  corrections  from  the  latent  heat,  hv,  differ  by  less  than  0.05  percent 
for  the  stated  variations  in  exterior  temperatures. 
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4  Coupled  Moisture,  Temperature,  and 
Structural  Analysis 


The  completion  of  the  hygrothermal  calculation  sets  the  stage  for  the  mechanical 
analyses  of  moisture  and  temperature  induced  stresses  in  the  FRP/concrete 
structure.  This  analysis  was  done  by  considering  the  moisture  and  temperature 
induced  strains  as  applied  forces  in  a  way  similar  to  the  application  of  body 
forces  as  input  parameters  in  standard  structural  analysis.  The  pre-defined 
structural  analysis  module  provided  in  the  FEMLAB  was  employed  to  carry  out 
this  coupled  hygrothermal-mechanical  analysis. 


In  the  coupled  hygrothermal-mechanical  problem,  the  total  apparent  strain  is 
the  sum  of  the  mechanical  strain,  the  temperature  induced  strain  due  to  thermal 
expansion,  and  the  moisture  induced  strain  due  to  moisture  swelling  effect.  The 
mechanical  strain  thus  can  be  calculated  by  the  following  equation: 


£  —  £ 

mechanical  total 


-£, 


thermal 


-£ 


moisture 


=  £ 


total 


-aAT-/3A9, 


Eq  9 


where: 

a  =  the  coefficient  of  thermal  expansion  (CTE) 

P  =  the  coefficient  of  moisture  swelling. 

The  constitutive  equation,  which  relates  the  stress  and  strain  states  in  a  solid 
medium,  can  then  be  expressed  as: 

li^mechanical  }=[D]{etotal-<xAT-pAG},  Eqio 

where: 

M=  {CT»  {J-  T,y  r,-  T,-  f  is  the  stress  vs*tor 

[fl]  is  the  stiffness  matrix 

{£ ^1 }  —  \zx  £y  £z  Yxy  Yyz  Y Xz  F  total  strain  vector  and  can  be  related  to  dis¬ 

placements  by  the  displacement-strain  relation  as- 
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Eq  11 


Using  the  same  finite  element  formulation  approach  as  shown  in  Section  3.1,  the 
resulting  finite  element  equation  for  structural  analysis  incorporating  moisture 
and  temperature  induced  strains  results  in  the  following  form: 

mM=(f.KK)+  thermal  {Pmoisture  }*  Eq  12 


where: 


{u}=  k  uyl  uzl  Ux2  uyl  uz2  •••  Uxr  uyr  uzrf  =  nodal  displacement 


vector 


[K]=  [\B]t[D\BW  =  element  stiffness  matrix 

Jv 

f  =  element  body  force  vector 

Jv 

{F5}=  f  element  force  due  to  surface  traction 

Js 

{^thermal }  =  I  l thermal)^ ^  =  element  force  vector  due  to  thermal  expansion 

Jv 

{Fmoisture}=\v\-BPyD^{Smoisture¥V  =  element  force  vector  due  to  moisture  swelling. 


In  the  above  expression,  [5]  =  {z}{TV}r  in  which  {Z,}  is  defined  in  Equation  6, 

? 

and: 


(/VO,y,z)}  = 


TV,  0  0  TV2 


0 


0 

0 


0  T^  0  0  TV2 

0  0  TVj  0  0  TV2 


TVr  0  0 

0  Nr  0 

0  0  TV 


Eq  13 


is  the  shape  function  matrix. 


Since  moisture  migration  in  the  structure  was  modeled  in  terms  of  humidity  po¬ 
tential  state  variable  instead  of  moisture  content,  the  swelling  coefficient,  p, 
which  is  generally  given  as  a  function  of  moisture  content,  needs  to  be  modified. 
The  moisture  induced  strain  in  terms  of  moisture  content  can  be  computed  as: 


Smoisture  =  OP  =  ~P, 


G_ 

P 


Eq  14 
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where: 


— 

0  [ — ]  =  the  specific  moisture  content  and  p  is  the  mass  density  of  the  material. 

kg 


The  3D  model  shown  in  Figure  4(a)  was  used  for  the  coupled  hygrothermal- 
mechanical  analysis  since  the  temperature  and  moisture  induced  stresses  in  the 
FRP/concrete  structure  need  to  be  described  in  a  three-dimensional  state.  A  fi¬ 
nite  element  mesh  as  shown  in  Figure  4(b)  was  generated  using  tetrahedral  lin¬ 
ear  elements.  The  density  of  the  mesh  was  chosen  based  on  the  considerations  of 
numerical  precision  and  the  required  computation  time  for  a  typical  nonlinear 
analysis  of  the  three-dimensional  model.  The  material  properties  of  the  concrete 
(Baker  1964)  and  FRP  (Tsai  1986)  are  given  in  Table  2,  where  the  swelling  coef¬ 
ficient,  |3,  is  approximated  as  a  constant  and  is  dimensionless,  representing 
strain  per  unit-specific  moisture  content. 

In  the  current  study,  a  nominal  dense  concrete  with  p  =  2140  kg/m3  and  a  typi¬ 
cal  glass/epoxy  FRP  with  density  p  =  1 160  kg/m}  were  used.  It  is  noted  that  the 
FRP  lamina  is  a  transversely  isotropic  material  and  the  numeral  subscripts  were 
used  to  indicate  its  three  principal  axes.  Along  the  direction  of  the  fiber  rein¬ 
forcement,  the  moisture-induced  swelling  is  negligible  since  the  glass  fibers  are 
not  permeable  and  have  a  much  higher  stiffness  than  the  epoxy  matrix. 

Figure  13  shows  the  schematic  of  the  model  with  prescribed  displacement  con¬ 
straint  boundary  conditions.  The  nodal  displacement  uz  on  the  all  nodes  on  the 

z  =  0  surface  is  constrained  based  on  the  consideration  of  symmetric  boundary. 
The  displacements  in  the  x  direction  for  the  nodes  on  the  central  edges  of  the 
model  with  y  =  -0.0194  and  y  =  0.194  are  also  restrained. 


Table  2.  Material  Properties  for  concrete  and  FRP. 


Concrete 

E 

(Gpa) 

V 

p 

a 

13.8 

0.18 

0.00767 

9.99  xIO-6 

FRP 

E2,  e3 

(Gpa) 

G12,G13 

(Gpa) 

G23 

(Gpa) 

V  V 

12,  13 

V23 

i 

ai 

a2,  a3 

38.6 

8.27 

4.14 

3.45 

0.26 

0.2 

0 

0.6  X  10'6 

8.6  x  10'6 

22.1  X  10'6 
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Figure  13.  Schematic  of  displacement  boundary  conditions. 

Finally,  to  prevent  rigid  body  motion  in  the  y  direction,  y-displacement  of  the 
node  at  the  coordinates  of  x  =  0.198,  y  =  0.194  and  z  =  0  is  also  restrained. 

The  indoor  and  outdoor  boundary  conditions  were  prescribed  as  Tin  =20 °C  and 
7 'mt  =  60°  C  for  temperature,  and  hin  =  50%  and  hout  =  97%  for  relative  humidity. 

Figures  14  and  15,  respectively,  show  the  resulting  steady-state  temperature 
and  relative  humidity  contours.  It  can  be  seen  that  the  humidity  and  tempera¬ 
ture  distributions  are  two-dimensional.  Also,  the  thermal  conductivity  of  the 
FRP  is  on  average  only  one  fifth  that  of  the  concrete,  a  higher  temperature  gra¬ 
dient  is  therefore  shown  in  the  FRP  laminate.  Because  the  FRP  laminate  only 
coves  the  central  portion  of  the  concrete  slab,  the  resulting  temperature  profile  is 
non-uniform  and  thus  the  heat  flux  flows  are  not  parallel.  It  can  be  seen  that 
the  maximum  temperature  gradient  in  the  FRP  region  exists  at  the  centerline. 
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Figure  14.  Temperature  contour  of  the  concrete  slab  partially  covered  with  an  FRP  laminate  (°C). 


Figure  15.  Relative  humidity  contour  of  the  concrete  slab  partially  covered  with  an  FRP 
laminate. 
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Both  the  humidity  and  thermal  humidity  diffusion  coefficients  for  concrete  are 
much  higher  than  those  of  the  FRP.  In  average,  they  are  2300  and  560  times 
greater,  respectively.  A  much  higher  humidity  gradient  therefore  is  shown  in 
the  FRP  region.  This  pattern  is  similar  to  that  of  the  temperature  distribution 
profile.  Again,  the  resulting  humidity  distribution  profile  is  non-uniform  and  the 
humidity  flux  flows  are  not  parallel.  The  FRP  laminate  behaves  as  a  moisture 
barrier  in  the  regards  of  moisture  diffusion. 

Figures  16  and  17,  respectively,  show  the  resulting  contours  for  the  interfacial 
nodal  shear  stress  T  for  the  concrete  and  FRP  regions.  This  stress  quantity 

represents  the  degree  of  shearing  mechanism  induced  on  the  surface  normal  to 
the  y-axis,  i.e.,  x-z  plane,  with  the  direction  along  the  x-axis.  It  is  the  stress 
component  responsible  for  interfacial  shear  debonding  in  the  transverse  direc¬ 
tion.  It  can  be  seen  that  the  extreme  value  of  this  stress  quantity  occurs  at  the 
discontinuous  interface  between  the  concrete  and  FRP,  and  that  the  stress 
gradually  returns  to  the  equilibrium  value  within  a  distance  that  is  five  to  six 
times  the  characteristic  FRP  thickness.  This  relaxation  phenomenon  is  a  gen¬ 
eral  feature  of  dissimilar  materials  in  intimate  contact,  subject  to  different  quan¬ 
tities  of  strain  and  is  known  as  “Saint-Venant’s”  principle. 

Figure  18  shows  the  contour  for  the  displacement  component  ux  on  the  deformed 
geometry  of  the  model  on  the  x-y  plane.  It  is  evident  that  the  concrete  slab  tends 
to  expand  more  than  the  FRP  due  to  the  combined  effect  of  thermal  expansion 
and  moisture  swelling.  As  a  result,  this  mismatch  in  the  expansion  in  the  con¬ 
crete  and  FRP  causes  the  interlaminar  shear  stress  in  the  x-direction.  Since  T 

shear  stress  acts  in  two  opposite  directions  on  the  two  edges  at  the  FRP-concrete 
interface,  the  stresses  along  these  edges  have  opposite  signs.  Some  shearing 
stresses  are  also  observed  on  the  side  surfaces  of  the  concrete,  i.e.,  at  the  x=0  and 
x=0.396  planes,  and  the  z=0.194  plane.  These  stresses  are  apparently  caused  by 
the  warping  deformation  about  the  z-  and  y-  axes,  respectively. 

Figures  19  and  20,  respectively,  show  the  resulting  contours  for  the  interfacial 
element  shear  stress  T  for  the  concrete  and  FRP  regions.  Based  on  the  avail¬ 
able  material  property  data,  the  FRP  has  no  moisture  swelling  effect  in  the  z  di¬ 
rection  (the  fibers  severely  constrain  swelling  in  this  direction),  while  the  con¬ 
crete  slab  expands  due  to  moisture  in  this  direction.  Mechanical  interaction 
mechanism  similar  to  that  discussed  above  also  occurs  in  the  z-direction,  result¬ 
ing  in  maximum  interlaminar  shearing  stress  T  at  the  interface  of  concrete 
and  FRP. 
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Figure  16.  Txy  shear  stress  contour  in  the  concrete  slab  (Pa). 


Figure  17.  Txy  shear  stress  contour  in  the  FRP  laminate  (Pa). 
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Figure  18.  Displacement  ux  contour  depicted  on  the  deformed  geometry  of  the  model  on  the  x-y 
plane  (m). 


Figure  19.  Tyz  shear  stress  contour  in  the  concrete  slab  (Pa). 
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Figure  20.  Tyz  shear  stress  contour  in  the  FRP  laminate  (Pa). 


In  addition,  the  nodal  displacement  uz  was  constrained  on  the  bottom  surface  to 
simulate  the  symmetric  boundary  condition.  No  interlaminar  shear  stress  Tyz  is 
induced  in  the  area  near  the  bottom  surface.  The  Tyz  stress  contour,  therefore,  is 
symmetric  with  a  mirror  image  about  the  bottom  surface. 

The  coupled  thermal  and  moisture  shear  stresses  rxy  and  vyz,  always  coexist  with 
other  mechanical  stresses  within  the  FRP-CMU  system.  Under  certain  circum¬ 
stances,  these  stresses  can  approach  the  tensile  strength  of  the  concrete  sub¬ 
strate.  A  common  engineering  rule  of  thumb  dictates  that  concrete  will  fail  in 
tension  at  about  10  to  15  percent  of  its  compressive  strength  (Ferguson  1979).  A 
typical  concrete  might  have  a  compressive  strength  on  the  order  of  4000  PSI 
(2.71  x  107  Pa),  and  thus  fail  in  tension  at  about  3.0  x  106  Pa.  Examination  of  zxy 
in  Figure  16,  indicates  that  the  maximum  concrete  stress  to  be  about  1.0  x  106 
Pa,  for  the  stated  extreme  conditions.  Laboratory  and  field  data  suggest  that, 
when  FRP  appliques  debond  from  a  CMU  substrate,  they  do  so  by  (tensile)  fail¬ 
ure  of  the  substrate.  Thus  the  finite  element  method  provides  a  valuable  metric 
for  the  qualification  of  what  particular  hygrothermal  circumstances  may  lead  to 
debonding  of  the  FRP  substrate. 

Finally,  since  all  the  results  presented  here  are  for  steady  state  conditions,  it  is 
helpful  to  get  some  measure  of  relevant  time  constants  associated  with  transient 
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changes  in  moisture  and  temperature.  Dimensional  analyses,  applied  to  equa¬ 
tions  1  and  3  readily  give  the  following  characteristic  time  constants  for  changes 
in  the  moisture  content  of  the  system: 

Te=d2IDe 

Here  d  and  De  are  the  material  thickness  and  moisture  diffusion  constant  (Ap¬ 
pendix  A),  respectively.  Using  characteristic  values  of  total  moisture,  0,  from 
Figure  9  for  CMU-FRP  and  d=  2  cm  for  CMU  and  0.2  cm  for  FRP,  one  finds  that 
Te  ~  100  hrs  for  CMU  and  Te  ~  730  hrs  for  FRP.  A  similar  analyses  from  Equa¬ 
tion  3  yields: 

^thm  ficdr / kf 

Again  using  characteristic  values  of  total  moisture,  Q,  from  Figure  9  for  CMU- 
FRP  and  d  -  2  cm  for  CMU  and  0.2  cm  for  FRP,  one  finds  that  Tthm  ~  5.2  min  for 
CMU  and  Tthm  ~  17.6  sec  for  FRP. 

With  these  characteristic  time  constants  for  the  moisture  and  thermal  change  of 
a  system,  it  is  possible  to  estimate  a  characteristic  strain  rate  of  change,  de/dr, 
for  a  given  unit  step  change  in  the  moisture,  0,  and  temperature  T.  Using  typical 
values  of  /?  and  a  from  Table  2,  one  finds  a  moisture  induced  strain  rate: 

ds/dz  ~  3.5  x  10'8  (strain)  /  hour  for  concrete,  and 
ds/dz  ~  5.13  x  10' 13  (strain)  /  hour  for  FRP. 

Similarly,  a  corresponding  thermal  induced  strain  rate  is: 

ds/dz  =  1 .2  x  10'4  (strain)  /  hour  for  concrete,  and 
ds/dz  ~  4.4  x  10'3  (strain)  /  hour  for  FRP. 

FRP  is  then  seen  to  contribute  the  dominant  transient  thermal  strain  to  the  sys¬ 
tem,  and  a  secondary  transient  contribution  derives  from  the  concrete  (CMU). 
These  strains  coexist  as  the  system  evolves  in  time  depending  on  prevailing 
boundary  conditions.  And  a  corollary  derived  from  these  dissimilar  rates  of 
strain  change  is  directly  applicable  to  experimental  design.  Any  simulation  at¬ 
tempting  to  emulate  natural  changes  in  the  environment,  must  take  into  account 
the  global  changes  in  the  CMU-FRP  system  that  will  have  vastly  differing  time 
constants,  depending  on  the  actual  physical  parameters  and  size  of  a  given  sys¬ 
tem.  This  further  underscores  the  need  for  an  explicit  coupled  hygrothermal  cal¬ 
culation  utilizing  the  finite  element  methods,  to  resolve  the  combined  strain  con¬ 
tribution  in  the  CMU-FRP  system. 
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5  Conclusions 


This  work  has  developed  a  generic  finite  element  formulation  for  coupled  hy- 
grothermal-mechanical  analysis  for  layered  structures.  The  theoretical  model 
follows  that  of  Philip  and  De  Vries,  for  mass  and  heat  transport  in  a  porous  me¬ 
dia  and  is  used  as  the  basis  for  the  present  work.  To  model  moisture  and  tem¬ 
perature  transport  in  structures  with  differing  physical  properties,  an  appropri¬ 
ate  continuous  moisture  representation  is  required.  This  is  because  the  total 
moisture  content  state  variable,  given  in  the  mass  and  heat  governing  equations 
is  discontinuous  at  the  interfacial  surfaces.  In  this  work,  the  relative  humidity 
of  the  surrounding  air,  which  is  related  by  the  sorption  isotherm  relations  to  to¬ 
tal  moisture,  was  chosen  as  the  continuous  variable  and  resulted  in  a  set  of  cou¬ 
pled  governing  equations. 

Based  on  these  modified  governing  equations,  a  finite  element  formulation  was 
developed  and  implemented  using  the  FEMLAB  finite  element  analysis  package. 
Solutions  to  this  set  of  coupled  equations  resulted  in  2D  level  curves,  represent¬ 
ing  the  thermal  and  moisture  distributions  within  the  FRP-CMU  system  plus  its 
defining  boundary  conditions.  The  resulting  moisture  and  temperature  induced 
stresses  in  the  FRP-CMU  structure  were  then  analyzed  using  the  humidity  and 
temperature  distributions  as  parameters  to  define  the  body  force  inputs  in  a 
structural  analysis  program.  To  accomplish  this,  the  predefined  structural 
analysis  mode,  provided  in  the  FEMLAB  software  package,  was  modified  and 
employed  to  carry  out  the  structural  analysis.  The  resulting  interfacial  shear 
stresses  caused  by  the  mismatch  of  the  thermal  expansion  and  moisture  swelling 
between  the  concrete  and  FRP  were  obtained. 

Continuing,  some  conclusions  derived  from  the  documented  research  literature 
and  simulations  conducted  in  this  research  effort  are: 

1.  When  FRP  is  applied  to  interior  CMU  surfaces,  it  is  known  from  laboratory  stud¬ 
ies  that  for  the  sake  of  both  economy  and  mechanical  response,  unidirectional 
FRP  reinforcement  in  the  form  of  fabric  strips  is  preferable  to  fabric  that  covers 
the  whole  surface  of  masonry  walls  (Triantafillou  1998;  Rivera,  Karbhari  2002), 
particularly  if  the  FRP  applique  is  deployed  in  extremes  of  humidity  and  tem¬ 
perature. 


Further,  this  type  of  applique  is  to  be  applied  to  clean,  unpainted  CMU  so  as  to 
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enhance  both  surface  adhesion  and  free  flowing  ventilation  of  the  CMU,  thus 
avoiding  moisture  pooling.  This  consideration  derives  further  support  from  the 
experience  of  exterior  paint  applications  in  cold  climes  (Korhonen  Bayer  1989), 
which  cautions  to  provide  ample  ventilation  of  painted  exterior  masonry,  to  avoid 
paint  delaminations,  primarily  from  efflorescence. 

Finally,  the  requirement  of  a  freely  ventilated  surface  can  be  particularly  crucial 
for  a  porous  system  containing  both  invading  and  receding  fluids  [Aker  1996; 
Descamps,  1997].  Water  imbibition  into  porous  building  materials  is  a  two-phase 
flow  process  in  which  the  imbibing  water  displaces  the  air  that  is  initially  present 
in  the  pore  space.  Accurate  non-destructive  measurements  of  transient  moisture 
profiles  during  imbibition  into  calcium  silicate  brick  indicate  that  the  air  outflow 
boundary  condition  has  a  predominant  influence  on  the  water  imbibition  process. 
As  a  consequence,  significant  differences  have  been  observed  in  the  water  sorp¬ 
tion  and  transport  coefficients  derived  from  water  imbibition  experiments  with 
free  and  bounded  air  outflow  boundary  conditions. 

2.  For  extreme  conditions  of  moisture  and  temperature  differentials,  such  as  en¬ 
countered  in  the  simulations  of  this  research  effort,  debonding  of  the  FRP  from 
the  CMU  substrate  is  possible.  This  comes  about  because  the  coupled  thermal 
and  moisture  shear  stresses  Txy,  always  coexist  with  other  mechanical  stresses 
within  the  FRP-CMU  system.  Under  certain  circumstances,  these  stresses  can 
approach  the  tensile  strength  of  the  concrete  substrate.  Laboratory  and  field 
data  suggest  that  when  FRP  appliques  debond  from  a  CMU  substrate,  that  they 
do  so  by  (tensile)  failure  of  the  substrate.  Furthermore,  any  experimental  or 
computer  simulation  attempting  to  emulate  natural  environmental  changes 
must  take  into  account  the  global  changes  in  a  CMU-FRP  system  that  will  have 
vastly  differing  time  constants,  depending  on  the  actual  physical  parameters  and 
size  of  a  given  system.  This  further  underscores  the  need  and  contributions  from 
an  explicit  coupled  hygrothermal  calculation,  using  finite  element  methods.  Do¬ 
ing  so  can  resolve  the  combined  strain  rates  in  a  given  CMU-FRP  system  that 
may  contribute  to  FRP-substrate  debonding. 

3.  The  boundary  conditions  used  in  the  mechanical  stress  simulations  of  this  report 
are  that:  The  inside  surface,  is  partially  covered  (60  percent)  with  FRP  and  is 
constrained  to  50  percent  RH  and  20  °C,  while  the  exterior  surface  is  held  at  97% 
RH  and  the  temperature  at  60  °C.  As  such,  the  results  derived  from  these  stress 
simulations  must  be  viewed  as  somewhat  extreme,  so  that  actual  Geld  conditions 
of  moisture  and  temperature  gradients  will  in  all  reality  be  less  severe. 

Given  this,  the  actual  simulated  shear  stresses  (item  2.  above)  were  less  than  the 
tensile  strength  of  the  underlying  CMU  by  at  least  a  factor  of  three.  This  conclu¬ 
sion  is  supportive  of,  and  consistent  with  a  decade  of  field  data  that  continues  to 
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justify  and  champion  the  use  of  FRP  upgrades  for  URM,  when  done  in  a  consis¬ 
tent  and  conscientious  manner  by  trained  professionals. 

4.  The  research  literature  cautions  that  the  true  kinetics  of  moisture  absorption 
(FRP/CMU)  under  freeze-thaw  exposure  in  aqueous  solutions,  is  not  well  under¬ 
stood  because  of  the  presence  of  microcrack  and  fiber-matrix  debond  related  deg¬ 
radation  and  possible  mass  dislocation  (Rivera  and  Karbhari  2002).  This  work 
considered  only  the  hygroscopic  transport  of  water  through  a  CMU-FRP  system 
for  temperatures  greater  than  0  DC  but  less  than  ~95  DC.  Given  this  limitation,  a 
natural  recommendation,  deriving  from  items  2.  and  3.  above,  is  to  extend  the 
present  modeling  and  simulation  effort  to  temperatures  below  the  freezing  point 
of  water.  In  this  temperature  range,  the  physical  properties  and  limitations  of 
epoxies  and  FRP  are  anticipated  to  dominate  those  of  the  CMU  substrate. 

5.  Appendix  A  to  this  report  fists  the  required  physical  parameters  that  appear  in 
the  defining  hygrothermal  Equation  set  1  to  3.  These  parameters  were  extracted, 
compiled  and  documented  from  the  available  (unclassified)  research  literature. 
Every  attempt  has  been  made  to  validate  and  cross  reference  these  parameters, 
to  ensure  their  correctness  and  validity.  This  validation  procedure  involved  the 
painstaking  checking  of  all  data  for  consistency  of  units  [MKS],  compatible  and 
relevant  usage.  In  some  cases,  differing  representations  of  the  same  physical 
constant  and  its  hygrothermal  representations  existed.  For  these  cases,  every  at¬ 
tempt  was  made  to  search  the  available  research  literature  and  seek  out  the  most 
consistent  and  well  founded  representation.  As  such,  the  Appendix  of  this  report 
represents  one  of  the  major  contributions  associated  with  this  work  unit,  and  can 
be  an  invaluable  resource  to  the  knowledgeable  practitioner  in  this  field. 
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Appendix  A:  FRP-CMU  Material 

Parameters 


Some  Material  Properties 

Equations  1  through  3  (pp  7,  8,  9)  outline  the  defining  physical  principles  govern¬ 
ing  moisture  migration  through  porous  structures.  However,  the  application  of 
these  principles  to  any  particular  system  requires  an  accurate  and  detailed 
knowledge  of  the  relevant  individual  physical  parameters  and  their  sensitivity  to 
moisture  and  temperature.  Generally,  such  data  has  not  existed  in  traditional 
handbook  format,  and  what  data  did  exist  was  incomplete,  unique  to  a  given 
task  or  environment,  and  non-uniform. 

Within  the  past  10  to  15  years,  the  building  physics  community  has  made  a  con¬ 
certed  effort  to  develop  experimental  techniques  to  evaluation  of  these  functions 
and  to  compare  predicted  results  with  available  field  data.  This  effort  has  been 
predominately  centered  in  Northern  Europe  and  Canada,  and  is  producing  re¬ 
sults  to  standardize  and  assemble  all  relevant  data  into  a  single  set  of  docu¬ 
ments.  The  culmination  of  this  effort  has  been  a  series  of  documents,  the  most 
recent  of  which  is  Annex  XXIV,  Task  3  (Kumaran  1996).  The  Annex  represents 
contributions  from  a  number  of  member  countries  and  reflects  the  diverse  needs, 
interests,  material  properties,  manufacturing  and  experimental  methodology  of 
individual  members.  Still,  the  Annex  represents  the  best  documented  summary 
of  unclassified,  common  building  material  hygro-thermal  properties  currently 
available  to  research  groups. 


Concrete 

A  brief  survey  of  the  various  concretes  listed  in  the  Annex  clearly  implies  that 
the  lighter  concretes  can  absorb  proportionally  more  moisture.  A  typical  CMU 
found  in  the  present  study  has  a  (dry)  mass  density  of  order  2140  kg/m3,  which  is 
characteristic  of  a  standard  concrete.  Drawing  upon  the  listed  Annex  data,  lit¬ 
erature  references  (van  Zijl  1999;  Janz  1997),  and  standard  engineering  practice 
it  is  possible  to  deduce  hygro-thermal  material  properties  appropriate  to  this 
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density  and  having  a  global  error  of  some  10  to  15  percent.  Equations  Al 
through  A3  express  these  data,  with  the  moisture,  9,  measured  in  units  of  kg/m3, 
relative  humidity  (RH)  expressed  a  fraction,  and  an  assumed  (unless  otherwise 
qualified)  temperature  of  20  °C: 

Concrete:  Isothermal  Absorption  (wetting)  Function  (Kumaran  1996) 

9  (kg/m3)  =  ax3  +  J bx2  +  cx  +  d  Eq  Al 

where: 

a' is  the  relative  humidity,  expressed  a  fraction,  and 
a  =  26.68 
b  =-247.75 
c  =  123.45 

d  =  0.1076,  for  2200  kg/m3  concrete. 

Concrete:  Isothermal,  Desorption  (drying)  Function  (Kumaran  1996) 

6  (kg/m3)  =  ax4  +  bx3  +  cx2  +  dx  +  e  Eq  A2 

where: 

vis  the  relative  humidity,  expressed  a  fraction  and, 
a  =  257.64 
b  =-426.95 
c  =202.88 
d  =77.823 

e  =  —0.0104,  for  2200  kg/m3  concrete. 

Concrete:  Wetting  Moisture  Diffusion  Function,  De 

Concrete  is  a  heterogeneous  mix  of  sand,  aggregate  (gravel),  water,  and  cement 
and  is  typically  characterized  by  its  mass  density.  Both  the  Annex  and  van  Zijl 
(1999)  list  comparable  values  for  Dq. 

Iog10  (Dq  [m2/s])  =  0.019426-10.5  (dimensionless)  Eq  A3 

for:  0.0  <  6  [kg/m3]  <  200 

Concrete:  Drying  Moisture  Diffusion  Function,  D0  (van  Zijl  1999) 
log10  (Dq  [m2/s])  =  f (9)  (dimensionless), 


f(0)  =  a96  +  b95  +  c04  +  d93  +  e92  +  f9  +  g 


Eq  A4 
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for:  0.0  <  G  [kg/m3]  ^  200 
where: 

a  =  1.0108823  E-12 
b  = -7.846242  E-10 
c  =  2.4277506  E-7 
d  = -3.7851524  E-5 
e  =  3.071331  E-3 
/  =-0.09861641 
g  =-8.9905065 

Concrete:  Thermal  Moisture  Diffusion  Function,  DT 

The  Annex  (Kumaran  1996)  lists  a  representation  for  the  Thermal  Moisture  Dif¬ 
fusion  Function,  Drp  However  a  more  common  and  useful  alternative  formula¬ 
tion  is  possible,  based  on  the  “WUFI-ORNL/IBP,  Hygrothermal  Design  Tool  for 

Architects  and  Engineers”  format  (Kuenzel  et  al.  2000).  In  this  representation 
the  thermal  flux  vector,  D^JT of  Equation  1  (p  7),  and  which  describes  the  mois¬ 
ture  vapor  migration  is  replaced  as  follows: 

DtVT  [kg/m2.s]  =  (5p(T>  I  p)V(hPs{Tf)  Eq  A5 


and 


PS(D  [Pa]  =  P0  exp(-(AH/R)[1/7- 1/373.15]) 


Here  p  is  the  (dimensionless)  vapor  resistivity,  and  is  a  physical  parameter  that 
provides  a  measure  of  moisture  vapor  transport  through  a  material  medium. 
The  vapor  resistivity  can  have  some  sensitivity  to  both  humidity  and  tempera¬ 
ture.  For  concrete  with  a  typical  density  of  order  2200  kg/m3,  the  Annex  (Kuma¬ 
ran  1996)  lists: 


1  Ip  = 
AH/R  = 

Po  = 

T 

h 


[0.0068  +  (8.21  x  10'5)exp(5.66h)  ] 

5201.08  [°K],  measured  for  liquid  water,  and 

101,325  [Pa],  the  pressure  of  a  standard  atmosphere,  also  equal  to  76  cm 
of  Mercury. 

the  absolute  thermodynamic  temperature,  measured  in  degrees  Kelvin, 
and 

the  relative  humidity  (RH),  expressed  as  a  fraction. 


The  Annex  lists  a  version  of  <5^(7),  however  an  equivalent  representation  is: 

5pHT)  =  1.83  x  10'10  (77  273.15  °K)0  81  ( PQ I P)  [kg/m.s.Pa]. 
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Here  Prepresents  the  ambient  atmospheric  pressure,  which  is  taken  to  be  iden¬ 
tical  to  PQ  for  this  modeling  effort.  SjiT)  is  then,  the  calculated  value  of  the  va¬ 
por  permeability  of  water  through  air,  evaluated  at  a  standard  atmosphere.  At 
this  time,  S^T)  is  not  known  to  have  been  experimentally  determined  for  air. 

In  this  representation  of  D^JT,  the  thermal  flux  vector  is  now  seen  to  have  both 
temperature  and  moisture  sensitivity,  through  the  APg(7)  function.  When  oper¬ 
ated  upon  by  the  V,  both  h  and  T cross  terms  are  generated  containing  the  gradi¬ 
ent  of  the  complementary  term.  The  Pg(T)  V  (h)  term  will  (or  can  be)  combined 

with  the  Z^V^term  of  Eq  1  (p  7),  resulting  in  a  weak  thermal  sensitivity. 

Concrete  Thermal  Data 

The  specific  heat  capacity  of  concrete,  c,  is  known  to  have  some  temperature  sen¬ 
sitivity,  that  is  of  order  3  to  4  percent  for  temperatures  between  10  and  37  °C, 
and  was  derived  for  fully  saturated  normal-weight  samples.  A  similar  sensitiv¬ 
ity  may  exist  for  unsaturated  samples,  however  no  data  appear  to  exist  to  sup¬ 
port  this  assumption  (Van  Geem  and  Fiorato  1983).  Given  this,  the  Annex  list 
the  thermal  data  as  expressed  in  Equation  A6. 

Concrete  Volumetric  Heat  Capacity 

pc  [J/m3.°C]  =  p0c0  +  4187  [J/kg.  °C]0  EqA6 


where: 

pQ  =  2200  [kg/m3],  dry  material 

cQ  =  940  [J/kg.  °C]  of  dry  material 

6  =  moisture  content  in  [kg/m3] 

Thermal  Conductivity  of  Moist  Concrete 
kT  [Watts/m.  °C]  =  2.74  +  0.00320 

Polymers-Epoxy  and  FRP 

Reliable  data  on  polymers,  with  and  without  fibers,  is  even  more  tenuous  than 
that  of  traditional  building  materials,  and  must  be  gleaned  by  hand  from  avail¬ 
able  journal  literature.  An  Epoxy  resin  being  considered  for  use  in  this  study  is 
EPON  828RS  (Shell  Chemicals),  a  liquid  difunctional  bisphenol 
A/epichlorohydrin-derived  resin,  of  average  molecular  weight  378  grams/mol. 
The  corresponding  FRP  is  a  Vinyl  Ester-glass  fiber  composite,  EXTREN®,  for 
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which  there  exists  moisture  uptake  as  a  function  of  ambient  relative  humidity. 
In  either  case,  the  relative  permeability  to  moisture  is  small,  roughly  that  of  a 
Type-1  (polymer)  vapor  barrier,  and  some  two  orders  of  magnitude  less  perme¬ 
able  than  a  typical  CMU. 

Epoxy  Resin:  Moisture  Data 

Some  moisture  data  has  been  determined  for  Epoxy  by  direct  immersion  in  sev¬ 
eral  penetrant  liquids  of  interest,  including  distilled  water  (Chin  et  al.  1999).  No 
humidity  data  is  available,  and  it  is  assumed  that  the  vapor  sensitivity  is  negli¬ 
gibly  small  in  keeping  with  the  previous  remarks.  From  the  data  of  (Chin  et  al. 
1999),  one  can  deduce  the  Arrhenius  temperature  dependence  and  activation  en¬ 
ergy,  A H,  for  both  (assumed)  wetting  and  drying  moisture  transport  as  follows: 

De  [m2/s]  =  D0  exp  (-AH/RT)  Eq  A7 

Note  that:  Dr[kg/m.s.  °K]  is  unknown,  and  assumed  =  0.0;  see,  FRP  below. 


Where: 

AH  =  69742.97  [J/mol]  for  Epoxy,  and 

Do  =  0.11854  [m2/s],  the  diffusion  coefficient  in  the  reference  state, 

T  =  the  absolute  thermodynamic  temperature,  measured  in  degrees  Kelvin. 

R  =  8.314  [J/mol.  °K],  the  thermodynamic  gas  constant. 

FRP  and  Vinyl  Ester  Resin:  Moisture  Data 

A  limited  amount  of  actual  moisture  diffusion  data  is  available  for  the  Vinyl  Es¬ 
ter-glass  fiber  composite,  EXTREN®  (Verghese  et  al.  1999).  Measurements  were 
performed  for  relative  humidities  (RH)  of  75,  85,  and  100  percent  water  bath  en¬ 
vironments,  with  an  assumed  temperature  of  60  °C.  Based  on  the  corresponding 
moisture  saturation  levels  and  an  assumed  matter  density  for  a  polymer  (epoxy, 
FRP)  of  1160  kg/m3,  Verghese  et  al.,  assign  the  following  functional  relation  be¬ 
tween  moisture  density  and  relative  humidity. 

Matrix  and  FRP:  Isothermal  Absorption  Function 

e  [kg/m3]  =  7.53  (x)43  Eq  A8 


Where: 

vis  the  relative  humidity  (RH),  expressed  a  fraction,  and 

the  drying  function  is  assumed  to  be  identical  to  the  wetting  function. 
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Matrix  and  FRP  Moisture  Diffusion  (wetting/drying)  Function ,  De 

Continuing  with  the  results  of  Verghese  et  al.  (1999)  and  Chin  et  al  (1999),  one 
can  deduce  (approximately)  the  Arrhenius  temperature  dependence  and  activa¬ 
tion  energy,  AH,  and  moisture  functionality  for  both  wetting  and  drying  moisture 
transport  as  follows: 

D0  [m2/s]  =(G(0))[  Dq  exp(-AH/RT)] 


Where: 

AH  =  21832.8  [J/mol]  for  Vinyl  Ester  (Chin  et  al.  1999),  and 

Dq  =  exp  ( +AH/R1 ),  for  T  =  (273  +  60)  °K  (Hence,  Dq  =  2659.7,  the  dimen¬ 

sionless  diffusion  coefficient  premultiplier  in  the  reference  state;  here 
Do  is  chosen  so  that  Dq\s  consistent  with  Verghese  et  al.  (1999),  for  T 

=  333  °K). 

T  =  the  absolute  thermodynamic  temperature,  measured  in  degrees  Kelvin. 


The  quality  of  the  moisture  sensitivity  data  is  subject  to  errors,  limitations  of 
graphical  transcription  and  measurement;  as  such,  a  simple  linear  fit  is  most 
appropriate  to  the  data  (Verghese  et  al.  1999). 

G(0)  =  (ad  +  b)x  10'13  (m2/s) 


for: 

0.0  <  Q  [kg/m3  of  moisture]  <  10 
Where: 

a  =  -1.865297  [m3/kg] 

b  =  49.124578 

Matrix  and  FRP:  Thermal  Moisture  Diffusion  Function,  DT 

Actual  FRP  data  for  the  thermal  moisture  diffusion  function,  D ^  does  not  ap¬ 
pear  to  exist  at  this  time.  However  adopting  the  same  format  as  for  concrete, 
based  on  the  “WUFI-ORNL/IBP,  Hygrothermal  Design  Tool  for  Architects  and 
Engineers”  formulation  (Kuenzel  et  al.  2000),  permits  a  bounding  approximation. 

DtVT  [kg/m2.s]  =  (5p(T)  I  \s)V(hPs(T))  Eq  A10 


and 


PS(T)  [Pa]  =  PQ  exp(-(AH/R)[1/T- 1/373.15]) 
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li  is  the  (dimensionless)  vapor  resistivity,  and  is  a  physical  parameter  that  pro¬ 
vides  a  measure  of  moisture  vapor  transport  through  a  material  medium.  The 
vapor  resistivity  can  have  some  sensitivity  to  both  humidity  and  temperature, 
and  is  of  order  40,000  for  Epoxy  resin  (Kunzel  1995).  Sp{T)  and  all  other  physi¬ 
cal  constants  and  parameters  are  as  listed  in  Eq  A5. 


Epoxy  Polymer  and  FRP  Thermal  Data 

The  thermal  conductivity  of  a  widely  used  thermoset  resin,  Hercules  3501-6,  has 
been  measured  over  a  wide  temperature  range  (Chern  et  al.  1993).  In  the  deri¬ 
vation,  these  authors  assume  that  the  volumetric  heat  capacity  of  the  epoxy  is 
comparable  to  that  of  a  heating  resistance  wire,  say  copper.  Assuming  this  and  a 
typical  matter  density  for  a  polymer  (epoxy,  FRP)  of  1160  kg/m3,  one  then  has 
the  following  thermal  data: 

Epoxy  Polymer  Volumetric  Heat  Capacity 

pEcE  [J/m3  °C]  =  pQc0  +  41 87[J/kg  °C]0 

where: 

p0cQ  =  0.36  x  107  [J/m3  °C],  dry  material 
9  =  moisture  content  in  [kg/m3]. 

Thermal  Conductivity  of  Host  (Epoxy  Polymer)  Material  (100%  cured  and  as¬ 
sumed  dry  [Chern  et  al.  1993]): 

kT [Watts/m  °K]  =  aT4  +  bT3  +  cT2  +  dT  +  e  Eq  All 


where: 

T  =  the  absolute  thermodynamic  temperature,  measured  in  degrees  Kelvin 
and, 

a  =  -9.522  E-ll 

b  =  1.111  E-7 

c  =  -4.057  E-5 

d  =  4.263  E-3 

0.382 


e 
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Matrix  and  FRP  Thermal  Data 

The  thermal  conductivity  and  the  volumetric  heat  capacity  for  the  matrix  and 
FRP  composite  are  generalized  functions  of  the  individual  epoxy  and  fiber  glass 
constituents,  the  details  of  how  the  individual  applique  is  put  together  and  ap¬ 
plied  to  the  substrate.  Given  this  and  the  general  impracticality  of  in  situ  meas¬ 
urements,  the  next  best  approach  is  to  approximate  these  thermal  parameters 
following  standard  mixture  techniques  (Hashin  1968;  Halpin  1984;  Springer 
1981). 

FRP  Volumetric  Heat  Capacity 

The  use  of  mixture  techniques  to  derive  the  volumetric  heat  capacity  for  the  Ma¬ 
trix-Fiber  combination  gleans  support  from  the  fact  that  this  thermal  parameter 
is  related  to  the  global  volume  properties  of  the  combined  materials,  in  a  simple 

additive  way,  that  reflects  the  total  thermal  energy  per  unit  volume.  Thus  the 
epoxy  matrix-FRP  volumetric  heat  capacity,  is: 

Pe-F  cE-F  [J/m3.  °C]  =  VfPyCy*  (1— Vfj  Pece  EcI  A12 


the  total  Epoxy  Matrix-FRP  volumetric  heat  capacity  [J/m3.°C] 
the  Epoxy  volumetric  heat  capacity  with  moisture  functionality  from 
above 

p cfG]~  0-209  x  107  [J/m3.°C],  the  dry  volumetric  heat  capacity  of  glass 
fibers. 

However,  if  carbon  fibers  are  used,  then 

p%px  ~  P(fc=  0-128  x  TO7  [J/m3.  °C],  the  dry  volumetric  heat  capacity  of  car¬ 

bon  fibers,  and  vyis  the  fiber  volume  fraction,  around  0.60,  for  both 
cases. 

Thermal  Conductivity  of  FRP  Material 

A  particular  complication  for  the  Epoxy- FRP  composite  is  the  anisotropic  nature 
of  the  thermal  conductivity,  which  arises  from  the  unidirectional  fibrous  nature 
of  the  FRP  within  the  epoxy  matrix.  However,  it  has  been  suggested  that  the 
transport  property  of  unidirectional  composites  along  the  longitudinal  direction 
may  be  approximated  using  the  standard  rule  of  mixtures  approach  (Hashin 
1968;  Halpin  1984;  Springer  1981),  as  follows. 


where: 

pe-fce-f 

PfFe 

PxPx 


kL  [Watts  /  m.°C]  =  vf  kf  +  (1-Vf)  kT 


Eq  A13 
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where: 

kj^  is  the  longitudinal  thermal  conductivity  [Watts  /  m.  °C] 

Vf  is  the  fiber  volume  fraction,  around  0.60 

kf  is  the  thermal  conductivity  of  the  fiber,  1.05  Watts  /  m.  °C  for  glass  and  80 

Watts  /  m.  °C  for  carbon  in  the  longitudinal  direction  and  12.5  Watts  /  m.  °C  for 
carbon  in  the  transverse  direction 

krp  is  the  thermal  conductivity  [W  /  m.  °C]  of  the  epoxy  matrix  from  Equation  All. 

The  transverse  thermal  conductivity  of  a  continuous  fiber  FRP,  k which  is 

also  the  out-of-plane  conductivity  value,  uses  the  Halpin-Tsai  equation  Agraal 
and  Broutman  (1990): 

kTran 1  kT  =  ^  +  W  /  ^  D  W  Ec1 A14 

where: 

1  =  HkflkT)-^l[(kflkT)  +  1]  EqA-15 

and 

k p  v>are  defined  from  above 

and  it  assumes  that  circular  cross-sections  exist  in  the  fibers. 
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